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ABSTRACT 


This  thesis  details  the  application  of  the  Porter  Method  to 
the  design  of  a  digital  tracker  controller  for  the  FPCC  aircraft. 
The  aircraft  is  simulated  with  software  models  supplied  by 
Lockheed,  and  incorporates  horizontal  and  vertical  canards,  along 
with  maneuver  flaps,  ailerons,  rudder,  and  jet  flaps  for  vectored 
thrust . 

Separate  lateral  and  longitudinal  designs  are  accomplished 
at  0.15  mach,  sea  level.  Longitudinal  designs  for  0.6  mach,  sea 
level,  and  0.9  mach  at  30,000  ft.  are  also  completed.  For  the 
one  lateral  design  the  following  maneuvers  are  performed: 

Yaw  pointing 
Flat  turn 

Horizontal  translation 

For  the  three  longitudinal  designs  the  following  maneuvers  are 
performed : 

Pitch  pointing 
Direct  lift 
Vertical  translation 

One  coupled  maneuver,  a  jink  maneuver  incorporating  both  vertical 
and  horizontal  translation  demonstrates  the  combined  lateral  and 
longitudinal  controller  at  0.15  mach.  A  "universal"  controller 
is  found  for  the  0.6  and  0.9  mach  flight  conditions  and  is  tested 
with  the  same  longitudinal  maneuvers.  All  of  the  controllers  are 
again  tested  with  the  addition  of  a  one  sampling  period  delay 


included  between  the  controller  and  the  control  input  to  the 
actuators.  This  represents  the  time  delay  that  may  be  introduced 
in  the  controller  by  the  microprocessor  implementation.  Finally, 
the  controllers  are  minimized  in  regards  to  numbers  of  gains 
required  in  the  control  law  matrices  by  finding  common  elements 
that  can  be  set  to  zero. 
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CHAPTER  I 
INTRODUCTION 


Background 

Digital  computers  and  electical  wiring  are  rapidly  replac¬ 
ing  analog  controllers  and  mechanical  linkages  in  the  design  of 
aircraft  flight  control  systems.  This  is  true  for  both  military 
and  civilian  aviation. 

The  digital  flight  control  system  must  handle  multiple 
inputs  from  the  pilot  and  also  command  multiple  outputs  simultan¬ 
eously,  ensuring  optimal  aircraft  performance  for  varied  flight 
conditions.  The  digital  flight  control  system  (DFCS)  can  control 
more  flight  control  surfaces  than  a  pilot,  and  can  also  control 
them  simultaneously.  Sensing  elements  such  as  the  gyros  and  the 
accelerometers  allow  the  DFCS  to  use  detailed  Information  about 
the  aircraft's  attitude  to  command  flight  control  surface 
deflections  more  efficiently  than  a  pilot. 

Using  a  DFCS  allows  the  pilot  to  command  a  particular 
maneuver,  instead  of  the  pilot  having  to  control  the  flight 
control  surfaces  directly.  Thus,  the  pilot's  workload  is 
reduced,  increasing  pilot  efficiency  and  allowing  the  pilot  more 
time  to  concentrate  on  other  vital  tasks.  In  the  combat  arena, 
with  the  corresponding  pressures  that  are  applied  to  the  pilot, 
the  possibility  of  error  should  be  reduced  and  the  chances  for 
survival  should  increase. 

The  DFCS  generates  the  signals  necessary  to  deflect  the 


ft 


control  surfaces  in 
pilot.  Consequently, 


i 

$ 


response  to  a  commanded  maneuver  from  the 
classical  single  input  single  output  (SISO) 
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design  techniques  are  not  well  suited  to  this  multiple  input 
multiple  output  (MIMO)  problem.  Other  techniques,  such  as  LQG, 
that  require  all  states  to  be  either  accessable  or  estimated  lead 
to  complicated  calculations  due  to  the  lack  of  full  state 
information  on  the  typical  aircraft.  An  alternative  design 
method  has  been  developed  by  Professor  Brain  Porter  and  associ* 
ates  at  the  University  of  Salford,  U.K. 

The  Porter  method  proposes  a  straight  forward  approach  to 
the  MIMO  control  .system  design  that  does  not  require  complicated 
mathematics  nor  does  it  require  considering  each  output-input 
pair  separately  (Ref  1).  Previous  AFIT  master's  theses  have 
demonstrated  the  use  of  the  method  (Ref  2  and  3). 

A  hypothetical  aircraft  called  the  FPCC,  for  flight  propul¬ 
sion  control  coupling,  has  been  proposed  by  Pratt-Whitney ,  Lock¬ 
heed,  and  Honeywell.  The  aircraft's  main  feature  is  the  use  of 
thrust  vectoring  using  jet  flaps,  hence  the  FPCC  label.  The  air¬ 
craft  also  Incorporates  horizontal  and  vertical  canards,  along 
with  maneuver  flaps.  The  FPCC  is  longitudinally  statically  un¬ 
stable,  adding  to  the  aircraft's  maneuverability  (Ref  4). 
Resulting  from  these  features  is  the  aircraft's  ability  to  per¬ 
form  decoupled  six  degree  of  freedom  (6D0F)  maneuvers. 

Design  of  a  DFCS  for  this  aircraft  is  aided  by  two  computer 
programs,  a  flight  simulation  program  for  the  FPCC  provided  by 
the  contractors  (Ref  5),  and  an  interactive  design  program  incor¬ 
porating  the  Porter  design  method  called  MULTI  (Ref  6). 


This  Chesls  extends  previous  design  work  done  on  a  DFCS  for 
the  FPCC  aircraft  (Ref  2),  by  developing  3  multi-variable  tracker 
control  laws  for  each  of  3  flight  conditions,  using  the  Porter 
method.  Each  controller  performs  coupled  or  decoupled,  lateral  or 
longitudinal  maneuvers.  Attempting  to  demonstrate  robustness,  a 
universal  controller  for  two  of  the  flight  conditions  is  desig¬ 
ned,  and  Its  performance  evaluated  compared  to  the  individual 
controllers . 

In  an  attempt  to  partially  validate  the  short  takeoff  and 
landing  (STOL)  capability  of  the  FPCC,  the  first  flight  condition 
is  at  0.15  mach,  with  a  low  angle  of  attack,  at  sea  level.  The 
other  two  flight  conditions  are  identical  to  those  from  Bausch- 
licher's  thesis,  0.6  mach  at  sea  level,  and  0.9  mach  at  30,000 
feet.  This  thesis  repeats  designs  at  these  flight  conditions 
because  the  longitudinal  aircraft  model  has  been  changed.  There¬ 
fore,  the  designs  for  these  flight  conditions  are  in  the  longitu¬ 
dinal  mode  only. 

The  responses  for  each  design  are  then  reaccomplished 
incorporating  a  time  delay  of  one  sampling  period  between  the 
output  of  the  controller  and  the  input  of  the  aircraft  plant 
model,  demonstrating  any  destabilizing  effects  from  this  lag. 
The  designs  are  redone,  as  necessary,  to  compensate  for  these 
effects.  The  simulations  are  also  repeated  with  certain 
controller  matrix  elements  set  to  zero,  demonstrating  possible 
simplifications  that  can  be  made  to  these  matrices. 

Additionally,  this  thesis  includes  two  improvements  to 
MULTI,  the  computer  program  used  for  controller  design  based  on 


Che  PorCer  method.  An  additional  option  was  added  that  calcu¬ 
lates  the  figures  of  merit  resulting  from  the  time  simulation 
used  in  the  program.  Specifically,  the  maximum,  minimum,  and 
final  values  are  displayed,  along  with  the  times  at  which  the 
maximum  and  minimum  occur.  The  option  then  allows  the  user  to 
use  the  default  value  of  within  2X  of  the  final  value,  or  the 
user  may  specify  the  value,  which  is  to  be  used  to  calculate  the 
settling  time. 

The  other  modification  to  MULTI  incorporates  the  option  of 
command  line  file  names  when  invoking  the  program.  The  files 
named  contain  system  and  design  data  to  be  used  within  the 
program. 


Current  Knowledge 

Papers  published  by  Professor  Porter  and  associates  detail 
the  design  technique  that  is  used  in  this  thesis  (Ref  6).  Also, 
previous  theses  have  both  used  and  explained  in  detail  the  Porter 
method.  This  thesis  is  the  first  to  investigate  the  STOL  capabi¬ 
lities  of  the  FPCC  aircraft,  and  extends  upon  the  work  done  by 
Jon  Bauschlicher  completed  in  December,  1982.  Bauschlicher's 
thesis  designed  multi-variable  controllers  for  flight  conditions 
reflecting  high  subsonic,  transonic,  and  supersonic  speeds.  The 
MULTI  modification  in  this  thesis  extends  the  work  started  by 
Douglas  Porter  as  published  in  his  thesis  of  December,  1981  (Ref 
3).  Bauschlicher's  thesis  and  the  reports  by  Lockheed  constitute 
the  background  Information  on  the  FPCC  aircraft  and  the  simula¬ 
tion  program  (Ref  4,5). 
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Computers  used  include  the  CDC  mainframes  (ASD)  at  Wright- 
Patterson  Air  Force  Base,  and  the  author's  microcomputer. 
Computer  support  facilities  used  include  the  RJE  (Remote  Job 
Entry)  sites  at  both  AFIT  and  the  Air  Force  Flight  Dynamics 
Laboratory  (AFFDL),  along  with  the  main  site  at  the  ASD  Computer 
Center. 


Overview 

Chapter  I  is  the  general  introduction,  with  Chapter  II  being 
an  introduction  to  the  Porter  method  of  controller  design. 
Chapter  III  is  a  description  of  the  FPCC  aircraft  with  a  summary 
of  the  aircraft's  behaviour  within  the  STOL  regime,  presenting 
previously  unpublished  tabular  analysis  of  the  aircraft's  re¬ 
sponse  to  flight  control  surfaces  at  slow  speed  and  low  altitude. 
Chapter  IV  describes  the  use  of  the  FPCC  simulation  program,  a- 
long  with  explanations  of  the  input  and  output  of  the  program. 

Chapter  IV  details  the  design  and  testing  of  the  controllers 
for  each  flight  condition.  Responses  of  the  controllers  to 
commanded  inputs  for  six  decoupled  and  one  coupled  maneuver  are 
given . 

Chapter  V  describes  the  design  of  a  universal  controller 
that  can  be  used  for  the  two  non-STOL  flight  conditions.  The 
effects  of  a  one  sampling  period  time  delay  in  the  output  of  the 
PI  controller  is  studied  for  all  the  flight  conditions  and  the 
universal  controller.  Finally,  certain  elements  of  the 
controller  matrices  for  the  STOL  flight  condition  are  set  to  zero 


to  study  the  possibility  of  reducing  the  number  of  gains  in  these 
matrices. 

Chapter  VI  summarizes  the  results  and  conclusions  of  this 
thesis,  including  insights  gained  from  the  use  of  the  Porter 
method  along  with  suggestions  for  its  future  use.  Chapter  VI 
also  details  the  assumptions  and  problems  associated  with  using 
the  FPCC  simulation  program,  and  gives  some  suggestions  for  its 
continued  use.  Also  in  Chapter  VI  are  ideas  for  possible  future 
additions  to  MULTI. 
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CHAPTER  II 

DESIGN  OF  LINEAR  MULTIVARIABLE  TRACKING  SYSTEMS  BASED  UPON 


SINGULAR  PERTURBATION  METHODS 


Introduction 


h>  Vv 


The  methods  used  In  this  thesis  for  the  design  of  the  con¬ 
trollers  were  developed  by  Professor  Brain  Porter,  University  of 
Salford,  England,  and  his  associates  (Ref  1,8,9).  The  references 
comprise  a  fairly  complete  documentation  of  the  development  of 
the  design  methodology  commonly  called  the  Porter  method,  but  two 
of  the  references  are  more  useful  as  summaries  of  the  method. 
The  first,  (Ref  1),  details  the  approach  for  both  regular  and 
irregular  known  plants  using  a  proportional  plus  integral  (PI) 
cascade  compensator  which  Is  either  analog  or  digital.  The 
second  (Ref  8),  briefly  explains  the  approach  for  an  unknown 
plant  with  an  analog  PI  compensator.  This  chapter  summarizes  the 
material  from  these  references.  Note  that  for  the  known  plants 
only  the  digital  controller  Is  considered,  since  that  Is  the  type 
of  controller  used  in  this  thesis  and  in  the  program  MULTI. 

If  additional  decoupling  and/or  reduction  in  the  initial 
undershooting  of  the  outputs  is  required,  then  a  technique  based 
upon  the  B*  approach  is  used,  which  is  also  included  in  this 
chapter.  If  the  B*  method  fails,  then  there  is  another  alterna¬ 
tive,  which  is  not  explained  in  this  thesis  but  is  referenced 
(Ref  9). 

Throughout  this  thesis  the  description  of  the  continuous 
time  plant  is  given  by: 

x( t)-Ax( t)+Bu( t)  (I) 


-vV 

v  V 


ft 


y ( t ) ■ C  x ( t)+Du( t) 

where 


A 

B 

C 

D 

and 


continuous  time 
continuous  time 
continuous  time 
continuous  time 


plant  matrix  (nxn) 
control  matrix  (nxm) 
output  matrix  (lxn) 
feed  forward  matrix  (lxm) 


(2) 


n  ■  number  of  plant  states 
1  -  number  of  outputs 

m  ■  number  of  inputs  (sometimes  called  p) 

The  PI  controller  performs  both  the  tracking  and  disturbance 

rejection  tasks,  and  as  stated,  is  the  discrete  time  type  for 

this  thesis.  The  equations  governing  the  controller  are: 

x[(k+l)TJ-  exp {AT}  *  x(kT)  +  jexp{AT}  *  Bdt  *  u(kT)  (3) 

o 

y [kT] ■  C  *  x(kT)  +  D  *  u(kT)  (4) 

where 


exp{AT}  ■  sampled  data  plant  matrix 
)exp{AT}  *  Bdt  -  sampled  data  control  matrix 
C  ■  sampled  data  output  marix 
D  ■  sampled  data  feed  forward  matrix 
T  ■  sampling  time 

Note  that  the  u(kT)  is  piecewise  constant  over  the  sample  period, 
resulting  in  its  removal  from  inside  the  integral.  The  separate 
names  designating  the  sampled  data  matrices  are  not  used  in  this 
thesis  since  these  marices  are  not  explicitly  used. 

Figure  1  shows  the  PI  controller  used  in  this  thesis.  Since 
the  computer  aided  design  program  MULTI  is  used  to  assist  in  the 
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design  of  the  controller,  the  placement  of  the  sampler  is  dicta¬ 
ted  by  MULTI's  implementation  of  the  controller.  MULTI  actually 
implements  two  samplers,  one  in  the  feedback  loop  from  the  plant, 
and  the  other  in  the  commanded  input.  Mathematically,  the  one 
sampler  shown  in  Figure  1  is  equivalent  to  the  other  two.  The 
zero  order  sample  and  hold  device  shown  is  indicative  of  the  way 
MULTI  passes  the  control  input  to  the  plant,  with  the  duration 
of  the  hold  controlled  by  an  option  within  MULTI.  The  integral 
action  shown  in  the  figure  is  currently  implemented  in  MULTI 
using  an  approximation  to  the  1/s  function. 

Since  the  choice  of  which  design  method  to  use,  Unknown, 
Known/Regular,  or  Known/Irregular,  is  based  upon  the  plant  para¬ 
meters,  the  rest  of  this  chapter  explains  how  the  choice  is  made 
and  summarizes  the  method  used  in  each  case.  At  the  end  of  the 
chapter  the  B*  method  is  explained. 


Unknown  Plant  (Ref  8) 

This  approach  applies  typically  to  an  industrial  process 
where  little  is  known  about  the  plant  model,  or  in  the  situation 
where  developing  the  state  model  is  not  desired.  The  problem  is 
to  develop  the  controller  without  knowledge  of  the  A,B,C,  or  D 
matrices. 

The  first  step  is  to  isolate  the  plant  and  determine  its 
steady  state  transfer  function  G(0).  This  is  possible  as  long  as 
the  plant  is  asymptotically  stable.  Since  the  transfer  function 


then 


G(>)-C(XI  -A)  B 

—  —  —  A  —  — 


G(0)— CA  B 


»vvv 


Since  the  addition  of  the  PI  controller  must  preserve  stabi¬ 


lity,  the  rank  of  G(0)  must  be  equal  to  1,  the  number  of  outputs. 
This  requirement  means  that  the  number  of  outputs  must  be  less 
than  or  equal  to  the  number  of  inputs,  and  that  G(X)  must  have  no 
transmission  zeros  either  at  the  origin  or  in  the  right  half 


plane . 


The  control  law  for  the  PI  controller  is: 


u(t)-«£  Ke(t)+eKz(t) 


wi  th 

e(t)«v( t)-y(t)  (7) 

z(t)-  |e(t)dt  (8) 

o 

where 

u(t)  -  control  input  to  the  plant 

**■  ■  ratio  of  error  to  the  Integral  of  error 

^  •  scaler  multiplier 

K  ■  the  controller  matrix 

e(t)  ■  error  vector 

v(t)  -  command  input  vector 

y(t)  ■  output  vector 

z(t)  ■  integral  of  error  vector 

Since  the  system  has  a  cascade  vector  integrator,  a  constant  step 
command  input  vector  yields  zero  steady  state  error  and  the 
output  vector  follows  the  input  vector,  hence  tracking  and  dis¬ 
turbance  rejection  result. 

Manipulation  of  the  above  equations,  along  with  the  fact  the 


steady  state  error  equals  zero  yields  the  following: 

K*G (0)  [G(0)G(0)T  f'  (9) 

which,  if  the  number  of  outputs  equals  the  number  of  inputs  means 

that  G(0)  is  square  and 

K-G(0 ) '  (10) 

The  above  development  can  be  extended  to  the  digital  PI  con¬ 

troller  by  writing  the  control  law  equations  as: 

u(kT)-T[Ko  e(kT)+K(  z(kT)  ]  (11) 

e(kT)*v(kT)-y(kT)  (12) 

z[ (k+1 )T]«z(kT)+Te(kT)  (13) 

where 

K  —K  (14) 

K-eK  (15) 


Note  that  Equation  (13)  is  an  approximation  to  Equation  (8). 

Following  the  same  basic  reasoning  up'-d  above  to  derive  E- 
quation  (9),  the  controller  matrices  K„  and  K,  can  be  found  from: 

K  -  {G(O)’  [G(0)G(0)T  f  ‘  }  (16) 
where  2  is  a  weighting  matrix  of  scalar  diagonal  entries  chosen 
by  the  designer. 

The  choice  of  the  diagonal  elements  of  £  and  of  the  sampling 
time  T  will  affect  the  output  response.  This  means  that  fine 
tuning  of  the  output  time  responses  can  be  accomplished  by  ad¬ 
justing  these  values. 

Unfortunately,  this  method  is  only  valid  for  a  plant  with 
negative  eigenvalues,  since  G(0)  must  be  stable.  This  precludes 
its  use  with  statically  unstable  aircraft  such  as  the  FPCC.  Ad¬ 
ditionally,  if  an  angle  is  used  in  the  state  vector  then  the  der¬ 
ivative  of  that  angle  (a  rate)  cannot  be  used  in  the  output 
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vector  since  that  would  generate  a  transmission  zero  at  the 
origin . 

Reference  2  discusses  the  effects  of  varying  the  parameters 
la  the  above  equations.  Increasing  an  element  in  the  £  matrix 
speeds  up  the  corresponding  output  response.  Multiplying  the  2E 
matrix  by  a  scalar  greater  than  one  increases  the  PI  controller 
Inputs  and  outputs,  while  multiplying  by  a  scalar  less  than  one 
decreases  the  inputs  and  outputs.  Decreasing  the  sampling  time 
T  increases  the  response  of  the  controller  outputs.  Changings 
varies  the  amount  of  over  or  undershoot  in  the  outputs. 

Known  Plant  (Ref  1) 

Figure  2  shows  the  block  diagram  of  the  plant  with  known  A, 
B,  C,  and  D  matrices.  Note  that  the  D  matrix  is  missing,  which 
is  due  to  a  requirement  that  the  system  representation  not  con¬ 
tain  a  D  matrix.  If  there  Is  a  D  matrix,  then  either  the  outputs 
or  state  representation  must  be  changed  to  eliminate  the  feed 
forward  matrix.  A  D  matrix  occurs  whenever  there  is  an 
acceleration  term  in  the  output  vector.  The  plant  equations  must 
be  expressed  in  the  following  format  to  use  the  Porter  method: 


with 


(18) 

The  A  matrix  of  Equation  (1)  has  been  partitioned  into  the  fol¬ 
lowing  sections: 


A  (n-l)x(n-l) 


Ata  (n-l)xl 

Aai  lx(n-l) 

A  lxl 
-22 


The  B  matrix  has  been  also  partitioned  as: 

B2  lxl  (must  have  full  rank) 
with  the  C  matrix  partitioned  as: 

C(  lx(n-l) 


Note  that  in  these  equations,  as  in  all  that  follow  (and  also  in 
MULTI),  the  number  of  outputs  equals  the  number  of  inputs,  i.e. 
l«m.  Therefore  the  representation  of  the  matrix  partitions  uses 
1  only. 

Equations  (1)  and  (2)  may  be  put  into  the  format  of  (17)  and 
(18)  with  a  similarity  transformation.  This  transformation 
yields  the  transmission  zeros  directly,  but  does  not  actually 
have  to  be  performed  to  apply  the  design  method. 

Implementation  of  the  Porter  method  requires  that  the  plant 
be  controllable  (A,B  be  a  controllable  pair),  and  also  observable 
(A,C  be  an  observable  pair),  along  with  the  requirement  that  the 
matrix 


have  full  rank  of  n+1  to  ensure  controllability  using  the  PI 
controller. 

After  these  requirements  have  been  satisfied,  the  next  step 
is  to  calculate  the  rank  of  the  first  Markov  parameter,  which  is 


m  m  »  m  w.i 


the  matrix  product  _CB  ("Ca  Ba) .  The  choice  of  which  design 
method  to  use,  either  regular  or  irregular,  is  based  upon  the 
rank  of  CB.  If  this  rank  is  less  than  1,  then  the  irregular 
method  must  be  used,  while  if  this  rank  is  equal  to  1  then  the 
regular  method  can  be  used. 


Known/Regular  Plant 

If  the  product  CB  (or  Ca§a)  has  full  rank,  then  the  plant  is 
considered  regular  and  C  B  can  be  inverted,  which  is  a  neces¬ 
sary  condition  to  use  this  method.  The  control  law  for  the 

digital  controller  is  given  by: 

u(kT)*(l/T)[K^  e(kT)+K(£<kT)]  (19) 

with 

e(kT)-v(kT)-y(kT)  (20) 

and 

z[(k+l)T]-z(kT)+Te(kT)  (21) 

The  controller  matrices  are  found: 

£c-ot<CaB;1j',2€  <22> 

K  «(C  BJ2€  (23) 

-i  -a-a  _  « 

where 

^  -  diagonal  matrix  {<r;  ,(£  ,  .  . .  ,<J^  } 
constrained  by 

-l<(l-ffj)<l 

Note  that  in  an  actual  design,  the  value  of  epsilon  can  be 

changed  to  obtain  a  required  range  of  values  for  the  elements  in 
the  controller  matrices,  and  thus  the  magnitudes  of  the  entries 


in  the  sigma  matrix  can  be  kept  within  the  required  range. 

Once  again  the  cascade  vector  Integrator  drives  the  steady 


state  error  vector  to  zero  for  a  constant  command  input  vector, 


eJV. 
■  V*V 


which  ensures  tracking  and  disturbance  rejection. 

As  the  sampling  time  T  approaches  zero  the  transfer  func¬ 
tion  matrix  G(A)  approaches  an  asymptotic  form: 

ra)-f<>)+ru>  c24> 


where 


?<*>-<*!„ -I. 


with 


(25) 

(26) 


GD 


A  -  -K  K 
—o  -a  -• 


A  C"'  if '  K  A  -A  c"  C 

—•a  -a  -  a  — «  —  n  —tx  - 2.  - 


2e-  I  5 


-I 

A  C 
.-•a  -a 


i] 


(27) 


(28) 

(29) 


Equation  (25)  corresponds  to  what  are  called  the  slow  modes  and 
(26)  corresponds  to  what  are  called  the  fast  modes.  Professor 
Porter  shows  that  as  T  approaches  zero  the  slow  modes  become  both 
uncontrollable  and  unobservable  while  the  fast  modes  remain  both 
observable  and  controllable.  Consequently,  as  T  approaches  zero, 
the  slow  modes  vanish  and  the  closed  loop  transfer  function 


becomes : 


r,(A)-ra)-(Ajfa-iJl-waBjKe)  cjjb2ko 

This  means  that  if  K0  is  chosen  such  that: 
CaBaKo-  diagonal  matrix  ,<T2  ,  .  . .  } 


then 


(30) 


(31) 


P(A)-diagonal  matrix  {07  /A-l+07,  (JZ/\-l  +  <ri . <7J  /A  -1  }  (32  ) 


V  *v’ _>*  /. /, /v,  •*  .  .  /  •*«*’-  *.  •*,  •*,  •*.  -  v-  % 


As  stated  earlier,  the  fast  mode  roots  and  the  transmission  zeros 


must  lie  within  the  unit  circle  to  ensure  stability.  The 
transmission  zeros  are  a  subset  of  the  slow  modes  and  are  given 


z4(X)-{|A  I  -I  -TA(t  +TA( 2  C*'  C(  |  -0} 
as  long  as  the  state  equations  are  in  the  form  of  (17). 


(33) 


From  Equation  (32)  the  conclusion  drawn  is  that  as  T  goes  to 
zero  decoupling  is  achieved.  However,  the  requirement  that  the 
transmission  zeros  lie  within  the  unit  circle  in  the  z  domain  may 
dictate  a  change  in  the  C  matrix.  This  means  a  redefining  of  the 
output  variables,  which  could  entail  a  major  change  in  the  ap¬ 
proach  to  the  design,  since  changing  C  must  not  affect  the  con¬ 
trollability  of  the  system. 


r 


Known/ Irregular  Plant  (Ref  1) 

When  the  product  CB  does  not  have  full  rank,  and  thus  Is  not 
Invertible,  the  plant  is  called  irregular.  In  this  case  extra 
plant  output  measurements  are  introduced  and  incorporated  with 
inner  loop  compensators.  This  results  in  the  feedback  vector: 


i(e)-jg,*At  c^aJ  f^(t> 

xa(t)j 


5,(0 


5a<fc> 


(34) 


(35) 


wl  th 


M  lx(n-l) 


F(  lx(n-l) 

F  lxl 
-z 


which  is  shown  in  Figure  3.  The  choice  of  the  measurement  matrix 
M  Is  such  that  FB  «  F^B^has  full  rank.  Basically,  the  design 
method  for  the  irregular  plant  is  just  like  that  for  the  regular 
plant,  with  F  replacing  C.  The  control  law  equation  for  the  PI 
controller  is  the  same  as  for  the  regular  plant,  but  the  error 
vector  becomes: 


e(kT)-v(kT)-w(kT) 


(36) 


Since 


then 


llm  [A„  5a(t)  ]■$ 

t-*oo  "* 


lim  e(kT)-lim  {v(kT)-y(kT)-M[A((  x  (  ( t )  +A|a  x^  (  t )  ]  } 
k-*°o  k  -*<*>  ”  '  '  * 

t-»oo 


-lim  { v( kT)-y ( kT) } 
k-»  oo 


Consequently  the  steady  state  error  vector  is  zero  for  a  con¬ 
stant  command  Input  vector.  Once  again  tracking  is  achieved. 


As  in  the  case  of  the  regular  plant,  the  closed  loop  trans¬ 
fer  function  matrix  approaches  an  asymptotic  value  as  the  gain 
1/T  goes  to  infinity,  of: 

P<X>-rOO+?<A>  (38) 


where 


?<>>-<£,  'Ml'  E,  ><*!.-«  l,  >" 


(T^,aEr>  ^39> 

(40) 


Once  again  the  r(X)  contains  the  slow  modes  and  P(\)  contains  the 
fast  modes.  The  transmission  zeros  are  given  by  (see  equation 


(33)): 

+TA,*/;'  P  |  -0)  (41) 

The  transmission  zeros  must  lie  in  the  left  half  side  of  the  s 

plane  or  within  the  unit  circle  in  the  z  domain. 

Unlike  the  regular  plant  case,  P  (X)  does  not  approach  zero 

as  the  sampling  time  decreases,  so  that  for  the  irregular  plant 
a 

both  TO)  andT  (X)  must  be  diagonal  for  decoupling.  The  choice 
of  the  measurement  matrix  M  and  the  control  law  matrices  Ko  and  K 
to  achieve  decoupling  requires  that  the  following  conditions  be 
satisfied : 

1.  M  is  chosen  such  that  Fa  and  F^B^  have  full  rank  •  1 

2.  All  closed  loop  poles  must  lie  in  the  left  half  plane 

3.  F  B  K  »(C„+MA.  )B  K  -diagonal  matrix  {<JT  ,<r.  ,  .  .  .  } 

4.  -1<(1-0T)<1  (roots  must  lie  within  unit  circle) 

The  standard  procedure  for  the  selection  of  the  measurement 
matrix  is  finding  the  sparsest  M  that  yields  an  F^  of  full  rank, 

-i 

while  also  diagonalizing  the  product  C  F  .  Although  this  tech- 


nique  is  straightforward,  it  is  not  always  possible  to  generate 
a  diagonal  transfer  function  matrix  (.()  . 

Another  problem  with  this  method  is  that  the  possible  pre¬ 
sence  in  the  output  of  the  modes  corresponding  to  the  transmis¬ 
sion  zeros  prevents  the  output  response  from  being  Improved 
beyond  a  certain  point  with  the  reduction  of  T. 

Once  the  measurement  matrix  is  chosen,  then  the  controller 
matrices  are  chosen  such  that: 


2* >  diaglcr;  ,<rA . ^  } 


(42) 


B*  Design  (Ref  1  and  9) 

This  design  technique  supplements  the  other  methods  when 
they  do  not  achieve  the  desired  decoupling  and/or  too  much  under¬ 
shoot  occurs  in  the  output  response.  Additionally,  under  cer¬ 
tain  conditions  the  B*  matrix  can  be  used  to  help  pick  the  ele¬ 
ments  in  the  M  matrix,  which  will  be  discussed  at  the  end  of  this 
section. 

The  first  step  is  to  form  the  B*  matrix  using  the  following 


formula: 


B*  - 


VaJ'  b 

Cr/L*  B 
— a  - 


(43) 


T  JjL 

c  A  B 


where  c  is  the  ith  row  of  C,  and  d;  is  the  smallest  integer  such 

—  t  — 

that  the  ith  row  of  B*  contains  at  least  one  non-zero  element. 
If,  however,  there  is  no  value  for  d^  that  yields  a  row  not  equal 


to  all  zeros,  then  d<  *n-l  and  the  ith  row  of  B*  is  set  to  all 

b  — “ 


zeros.  When  this  happens  B*  is  not  of  full  rank  and  this  method 
cannot  be  used.  However,  under  certain  constraints  another  ap¬ 
proach  can  be  used  and  is  discussed  in  the  last  paragraph  of  this 
section . 

When  B*  has  full  rank  the  controller  matrices  can  be  deter¬ 


mined  from: 


K  -B*  Z 

—  0  -  —  I 

K(-G(oj' 


(44) 

(45) 


where  5  and  2^  are  diagonal  lxl  matrices  chosen  by  the  designer, 
based  on  the  concepts  that  2  controls  the  initial  output  re¬ 
sponse  and  affects  the  steady-state  response.  As  with  the 


other  methods  the  control  law  equations  are: 
u(kT)-(l/T)  [K  m(kT)+K  s(kT)] 


(46) 


e(kT)-v(kT)-y(kT)  (47) 

If  B*  is  singular,  the  above  method  cannot  be  used,  but 
there  is  another  possibility.  If  the  row  dimension  of  the  C  ma¬ 


trix  (Equation  (18))  Is  greater  than  the  row  dimension  of  A^  (E- 
quation  (17)),  then  B*  can  be  used  to  help  pick  the  elements  in 
the  M  matrix  (using  the  Irregular  plant  method).  When  using  this 
approach,  the  B*  matrix  must  be  formed  by  using 

£-c, 

A-A(< 

in  Equation  (43). 

The  technique  is  continued  by  forming: 


F  -C  +MA,, 

— 2  —  •* 

and  assigning  values  to  those  entries  in  F  that  correspond  to 

*  a 


the  non-zero  elements  of  B  *  .  The  design  is  continued  along  the 
lines  of  the  irregular  plant  method,  with  further  refinement  of  M 
based  upon  the  fact  that  must  be  diagonal,  F  must  have  full 
rank  and  M  should  be  as  sparse  as  possible. 


Summary 


The  design  methods  presented  in  this  chapter  form  the  basis 
for  the  Porter  method,  as  developed  by  Professor  Brian  Porter  and 
his  associates.  The  computer  program  MULTI  incorporates  these 
concepts  into  an  interactive  tool  used  to  design  controllers  such 
as  those  designed  for  this  thesis.  Although  the  presentation  in 
this  chapter  specified  digital  PI  controllers,  since  that  is  how 
MULTI  i s  written.  Professor  Porter's  publications  include  similar 
design  techniques  for  analog  controllers. 


CHAPTER  III 


FPCC  AIRCRAFT/SIMULATION  PROGRAM  DESCRIPTION 

Introduction 

This  thesis  designs  controllers  for  the  FPCC  aircraft,  an 
alrcaft  which  exists  in  concept  only.  This  “paper  airplane"  was 
conceived  through  the  combined  efforts  of  Lockheed,  Pratt  and 
Whitney,  and  Honeywell.  FPCC  designates  Flight  Propulsion  Con¬ 
trol  Coupling,  which  summarizes  the  most  striking  feature  of  this 
plane,  the  ability  to  vector  thrust  with  the  jet  flaps.  This 
capability,  taken  with  the  other  flight  control  surfaces. 
Increases  the  control  vector's  dimensions,  allowing  for  more 
commanded  quantities  in  the  input  vector. 

The  contractors  developed  the  model  and  simulation  program 
between  June  1976  and  June  1977.  At  that  time  the  techniques  de¬ 
veloped  by  Professor  Porter  had  not  yet  been  published,  ana  the 
Air  Force  was  not  involved  in  designing  aircraft  flight  control 
systems  for  aircraft  capable  of  decoupled  six  degree  of  freedom 
maneuvers.  But  the  FPCC  aircraft  was  designed  with  the  extra 
flight  control  surfaces  needed  for  such  maneuvers,  so  that  as  the 
Air  Force  developed  an  Interest  in  decoupled  maneuvers  the  FPCC 
became  a  logical  plane  to  test  the  ideas  on.  In  keeping  with  the 
concept  of  using  the  FPCC  to  test  the  Porter  method  of  controller 
design,  this  thesis  extends  the  work  of  Jon  Bauschlicher  (Ref  2), 
continuing  to  use  the  maneuvers  similar  to  those  currently  being 
tested  on  the  AFTI  F-16  aircraft.  Additional  information  on  the 
FPCC  aircraft  can  be  found  in  References  4  and  10. 
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Aircraft  Description 

Figure  4  shows  the  FPCC  aircraft.  This  aircraft  is  a  single 
seat  supersonic  fighter  with  a  primary  mission  of  air  superiority 
fighter,  and  a  secondary  mission  of  ground  attack.  The  airframe 
has  the  following  dimensions: 

Aircraft  Length  -  60.0  ft 
Wing  Span  (B)  ■  54.25  ft 
Wing  Mean  Aerodynamic  Chord  ■  14.00  ft 
Wing  Planform  Area  (S)  *  654.00  ft 
Aspect  Ratio  (B*B/S)  *  4.50 
Takeoff  Mass  ■  1055.9  slugs 
Maneuver  Mass  *  900.62  slugs 
Maximum  Angle  of  Attack  -  +23  degrees 
Minimum  Angle  of  Attack  ■  -11.5  degrees 
The  FPCC  has  standard  flight  control  surfaces  in  aileron  and 
rudder,  along  with  non-standard  surfaces  in  its  horizontal  and 
vertical  canards  and  maneuver  and  jet  flaps.  The  vertical  ca¬ 
nards  along  with  the  rudder  give  the  FPCC  the  opposing  lateral 
moments  which,  when  combined  with  the  aileron,  allow  decoupled 
lateral  maneuvers.  Likewise,  the  jet  flaps,  maneuver  flaps,  and 
horizontal  canards  supply  the  opposing  longitudinal  moments 
necessary  to  perform  decoupled  longitudinal  maneuvers.  There  are 
no  elevators. 

Table  I  shows  a  summary  of  the  FPCC  aircraft  design  guide¬ 
lines  and  constraints,  and  Table  II  shows  the  control  surface 
limits  along  with  the  actuator  and  sensor  dynamics  for  the  con¬ 
trol  surfaces.  Note  that  although  these  are  the  same  actuator 
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SUMMARY  LISTING  0?  INITIAL  FPCC  REFERENCE  AIRCRAFT 
POINT  DESIGN  GUIDELINES  AND  CONSTRAINTS  (REF  a) 


1.  Primary  aircraft  design  mission  —  supersonic/transonic 
air  superiority 

Secondary  aircraft  design  mission  —  transonic  close  air 
support 

Aircraft  dash  speed  capability  —  within  the  range  of 
Mach  2.2  through  3*0 

2.  Advanced  CTOL  or  STOL  design  emphasizing  air  combat  maneu¬ 
vering  and  air  combat  tracking  requirements  for  high  sub¬ 
sonic/transonic  air-to-air  and  air-to-ground  tasks 

3.  Mission  range,  payload,  takeoff  and  landing  distance  re¬ 
quirements  are  to  be  considered  of  secondary  importance 

4.  Aircraft  and  propulsion  design  not  constrained  to  any 
specification  except  for  those  noted  in  the  RF?  Statement 
of  Work 

5.  Flight  control  designs  may  include  any  advanced  feature 
such  as  Relaxed  Static  Stability  (RSS)  and  Direct  Side 
Force  Control  (DSFC) 

6.  Aircraft  design  may  include  modulatory  type  Aerodynamic/ 
Propulsive  Interactive  Force  (APIF)  systems 

7.  Digital  control  and  fly-by-wire  assumed 

8.  The  supersonic  engine  air  inlet  shall  be  of  the  external  or 
mixed  compression  type  consistent  with  the  selected  air¬ 
craft  dash  speed  Mach  number.  Inlet  bypass  air  may  be  used 
for  aircraft  maneuvering  augmentation 

9.  Aircraft  design  thrust  to  weight  (T/W)  ratio  shall  be 
higher  than  1.0  with  dry  engine  operation.  However,  engines 
shall  Incorporate  afterburning.  Engine  control  may  Include 

-  variable  fan  guide  vane  angle 

-  variable  compressor  stator  angle 

-  variable  turbine  area 

-  variable  exhaust  nozzle  area 

Engine  bleed  capabilities  may  be  considered  consistent  with 
APIF  system 

10.  The  Lockheed  reference  aircraft  point  design  will  incor¬ 
porate  the  current  state  of  the  art  Pratt  and  Whitney  air¬ 
craft  F-100-PW-100  afterburning  turbofan  engine  and  controls 


dynamics  used  in  the  program  MULTI,  these  are  not  the  same  sensor 


dynamics.  MULTI  sensors  refer  to  the  sensors  used  to  sense  the 
aircraft  states,  not  the  flight  control  surfaces.  The  only  pos¬ 
sible  overlap  for  this  aircraft  would  be  if  the  total  thrust  were 
used  as  both  a  flight  control  and  a  state. 

The  original  design  for  the  aircraft  stipulated  that  the  jet 
flaps  would  be  limited  to  only  positive  deflection  (trailing  edge 
down)  of  0  to  90  degrees,  and  also  restricted  the  maneuver  flaps 
to  +  15  degrees.  Jon  Bauschlicher  had  the  simulation  program 

changed  to  reflect  the  values  listed  for  these  control  surfaces 
in  Table  II,  and  this  author  has  let  this  modification  stand. 

More  detailed  Information  about  the  workings  of  the  control 
surfaces  is  found  in  the  next  section  on  the  simulation  program. 

FPCC  Aircraft  Simulation  Program 

In  effect,  the  simulation  program  is  the  aircraft.  All  of 
the  control  surface  limits  depicted  in  Table  II  are  set  within 
the  program,  as  are  the  actual  implementations  of  these  surfaces. 
The  current  version  of  the  program  deflects  the  ailerons  and  rud¬ 
der  in  the  standard  manner  (ailerons  together  in  opposite  direc¬ 
tions),  and  deflects  the  maneuver  flaps,  the  horizontal  canards, 
the  vertical  canards,  and  the  jet  flaps  symmetrically.  The 
program  allows  asymmetric  thrust,  but  this  feature  is  not  used  in 
this  thesis. 

All  of  the  matrices  used  in  this  thesis  are  generated  by  the 
simulation  program,  as  are  the  initial  state  vectors.  The  simu¬ 
lation  program  is  available  from  the  Flight  Dynamics  Lab  at 
W r igh t-Pa t terson  Air  Force  Base,  and  was  developed  by  Lockheed 


TABLE  II 


CONTROL  SURFACE  LIMITS  AND  DYNAMICS 


Control 

Surface 

Deflection 
M<1. 0 

Limits 
M>1. 0 

Servo 

Dynamics 

Sensor 

Dynamics 

Maneuver  Flap 

<‘mp> 

+30° 

+15° 

25 

s+7* 

100 

s+TSff 

Jet  Flap 

±90° 

+90  • 

25 

100 

<V 

s+2!T 

s+loO 

Horizontal 

+20° 

+10° 

25 

100 

Canard 

s+25 

S+100 

Thrust 

0.574 

100 

*-3 

O 

•-3 

s+100 

Aileron 

+20° 

+20° 

25 

100 

<*a> 

(per  side) 

(per 

side) 

s+25 

s+100 

Rudder 

+30° 

+15° 

25 

100 

«r> 

s+25 

s+100 

Vertical 

+35° 

+15° 

25 

100 

Canard 

(6VC^ 

n 

s+25 

s+100 

*n  *>  *-9-  .5s  ■»  o*  *s  »3  »<D  ’  J2 


under  contract  to  the  Air  Force.  Running  in  batch  mode 
program  outputs  the  data  used  for  this  thesis  in  the  fo 
Equations  (48)  and  (49). 

p  -.696  .364  -.00818  0  000  0  0  0 

.0654  -.212  .00415  0  000  0  0  0 

31.5  -164  -.172  31.6  000  0  0  0 


-31.3  167  0  0 


0  -.103  .0134  -.000156  0 


0  166  -.638  -.215 


0  -31.5  .0867  -.0214 


-.982  .187 


.187  .982 


.00000326  .00000325 


-.000410  -.000415 


.0000443  .00000935 


0.0122 

0 

0 

.0118 

-.00212 

.00939 

. 0000322 

-.0000322 

0.00114 

0 

0 

-.0112 

.00220 

-.000238 

. 0000288 

-.0000288 

0 

0 

0 

.111 

.0363 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.0390 

0 

0 

0 

-.00908 

-.0000298 

-.0000298 

0 

-.111 

0 

0 

0 

-.118 

-.000367 

-.000367 

0 

-.0765 

0 

0 

0 

-.0211 

.00105 

.00105 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

The  output  equation  is 


0  0  0  0 


.00585  -.00112  0  0 


.00598  0000 


0  0  0  0 


.187  .982 


ACGX 


000  -.154  .0867  -.0214  0  .0000443 


ACGY 


.125  1.02  -.172  0000 


ACGZ 


000  .937  -.638  -.215  0  -.000410 


.00000935 


-.000415 


0  -.0765  0  0 


-.0211  .00105  .00105  -24.2  -24.2  6 


0  .111  .0363  0 


0  -.111  0  0 


-.118  -.000367  -.000367  0 


Note  that  these  equations  contain  the  output  for  the  trim 
condition  of  0.15  mach  at  sea  level,  which  is  the  first  flight 


condition  used.  Tables  III  and  IV  show  the  input  and  state 
vectors,  which  are  basically  the  same  as  those  used  in  this 
thesis,  while  Table  V  list  the  output  vector  which  is  not  the 
sane  as  used  in  this  thesis.  While  the  simulation  program 
generates  data  that  is  used  for  the  controller  design,  the  form 
of  the  matrix  equations  is  changed  slightly. 

The  simulation  program  actually  trims  the  aircraft  given  a 
mach  number,  altitude,  and  settings  for  the  maneuver  flaps  and 
jet  flaps.  The  plane  is  then  trimmed  using  the  horizontal  ca¬ 
nards  and  the  thrust.  The  ailerons,  vertical  canards,  and  the 
rudder  are  not  deflected.  Trimming  is  successful  only  if  the 
limits  on  the  flight  control  surfaces,  thrust,  and  angle  of  at- 


TABLE  III 

INPUT  VECTOR,  u« 


Variable 

Units 

Description 

6a 

deg 

sum  of  both  aileron  deflections 

6C 

deg 

horizontal  canard  deflection 

6mf 

deg 

maneuver  flap  deflection 

6r 

deg 

rudder  deflection 

5  vc 

deg 

vertical  canard  deflection 

6j 

deg 

Jet  flap  deflection  WRT  fuselage 
ref  line  ^ 

pi 

lbs 

Net  thrust  engine  1 

P2 

lbs 

Net  thrust  engine  2 

CDI1 

- 

Inlet  drag  coefficient  engine  1 

CDI2 

- 

Inlet  drag  coefficient  engine  2 

TABLE  IV 


STATE  VECTOR,  X* 


Variable 


Units 

Description 

rad/sec 

Body  axis  roll  rate 

rad/sec 

Body  axis  yaw  rate 

ft/sec 

Body  axis  side  velocity 

rad 

Roll  angle 

rad 

Yaw  angle 

ft 

Cross-range  position  re¬ 
ference  to  initial  body  axis 

rad/sec 

Body  axis  pitch  rate 

ft/sec 

Body  axis  vertical  velocity 

ft/sec 

Body  axis  forward  velocity 

rad 

Pitch  angle 

ft 

Altitude 

ft 

Down-range  position  re¬ 
ferenced  to  initial  body  axis 

TABLE  V 


OUTPUT  VECTOR,  y 


Variable 

Units 

Description 

a- 

rad 

Angle-of-attack 

6 

rad 

Sideslip  angle 

VEL 

ft/sec 

Total  airspeed 

ACGX 

ft/sec2 

Accel,  of  c.g.  along  x-body 
axis 

ACGY 

ft/sec2 

Accel,  of  c.g.  along  y-body 
axis 

Accel,  of  c.g.  along  z-body 
axis 

ACGZ 

ft/sec2 

tack  are  not  exceeded.  The  flight  conditions  used  in  this  thesis 


are : 


1). 


0.15  mach,  0  ft  altitude 


2) .  0.60  mach,  0  ft  altitude 

3) .  0.90  mach,  30,000  ft  altitude 

Since  the  first  condition  represents  a  landing  configuration,  the 
angle  of  attack  is  restricted  for  pilot  visibility  of  the  runway. 
The  simulation  program  trimmed  the  aircraft  to  11  degrees  angle 
of  attack  with  the  maneuver  flaps  set  for  +25  degrees  and  the  jet 
flaps  set  at  +15  degrees,  and  the  horizontal  canard  trimmed  to 
+0.5  degrees. 

The  second  flight  condition  represents  a  possible  ground 
attack  mode,  high  subsonic  speed  at  low  altitude.  The  third 
flight  condition  is  a  candidate  for  aerial  combat,  just  below 
supersonic  at  high  altitude. 

Table  VI  shows  the  results  of  several  different  flight 
conditions  run  on  the  FPCC  simulation  program.  These  conditions 
were  tested  to  observe  the  interaction  between  the  trimmed  angle 
of  attack,  horizontal  canard  setting,  thrust,  and  the  given 
settings  of  maneuver  flaps  and  jet  flaps  at  low  speeds  at  sea 
level . 

Table  VI  clearly  shows  the  relationship  between  the  angle  of 
attack  and  the  jet  flap  setting  as  slow  speeds.  However,  to 
counter  the  moment  caused  by  a  large  jet  flap  setting,  the 
horizontal  canards  must  be  deflected  possibly  beyond  their 
limits.  To  offset  this  drawback,  the  horizontal  canards  need 
more  effectiveness,  perhaps  through  increased  area.  This  would 
allow  for  a  smaller  angle  of  attack  at  even  slower  airspeeds  than 
used  in  this  thesis,  thus  permitting  a  more  STOL-like  operation. 

The  values  in  Table  VI  often  exceed  the  limits  listed  in 
Table  II  because  the  simulation  program  does  not  actually  impose 


TABLE  VI 

FPCC  AIRCRAFT  RESPONSE  AT  SEA  LEVEL 


Mach  # 

cX 

<fc 

£ 

F 

rror 

0.10* 

65.62 

-17.18 

30 

0 

16,668 

0.10* 

57.81 

-  1.79 

30 

10 

15,564 

0.10* 

51.42 

5.36 

30 

15 

15,366 

0.10* 

44.22 

12.34 

30 

20 

15,296 

0.10 

33.04 

27.16 

0 

30 

16,586 

0.10* 

30.17 

25.49 

30 

30 

15,210 

0.10* 

10.39 

74.00 

30 

60 

16,504 

0.15* 

13.44 

-  9.20 

30 

0 

8,826 

0.15 

13.44 

-  9.20 

15 

0 

8,826 

0.15* 

11.48 

-  2.52 

30 

10 

7,730 

0.15* 

10.78 

0.51 

30 

15 

7,468 

0.15* 

10.78 

0.51 

25 

15 

7,468 

0.20 

5.92 

-  2.31 

15 

0 

6,846 

i 

0.20* 

5.92 

-  2.31 

30 

0 

6,846 

0.20* 

4.84 

2.07 

30 

10 

6,702 

0.20* 

4.30 

4.26 

30 

15 

6,884 

*  Represents  a  flight  condition  with  a  column  of  all  zeros  in 
the  B  matrix  corresponding  to  the  maneuver  flaps. 


these  limits  on  the  aircraft  during  trimming.  However,  the 
program  indicates  whenever  the  limits  have  been  exceeded,  and 


then  continues  with  the  generation  of  the  linear  dynamics 


matrices  . 


Several  of  the  flight  conditions  yield  B  matrices  that  have 
a  column  of  zeros  in  the  column  corresponding  to  the  manuever 
flaps.  If  the  program  is  correct  in  generating  these  zeros,  then 
the  maneuver  flaps  have  no  effect  (for  small  pertubations )  at 
these  flight  conditions.  This  is  the  case  for  the  flight 
condition  at  0.15  mach  used  in  this  thesis,  which  is  why  the 
maneuver  flaps  are  removed  from  the  longitudinal  model  (Chapter 
IV).  However,  since  the  simulation  program  seems  to  be  erroneous 
in  certain  other  entries  within  the  B  matrix  (Appendix  A),  then 
suspicion  has  to  fall  on  all  of  the  B  matrix,  including  these 
columns  of  zeros.  The  position  taken  for  this  thesis  is  that 
these  columns  of  zeros  are  correct,  and  the  longitudinal  model 
adjusted  accordingly.  This  position  resulted  more  from  a  lack  of 
time  to  prove  otherwise,  than  any  other  consideration. 

An  interesting  aside  is  that  STOL  operation  at  very  slow 
speeds  is  going  to  require  large  amounts  of  thrust,  which  may  not 
be  a  desirable  trait. 

This  thesis  assumes  that  the  maximum  combined  thrust  is 
30,000  pounds,  based  upon  the  specification  of  at  least  a  1.0 
thrust  to  weight  ratio  and  the  weight  of  the  aircraft. 
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CHAPTER  IV 


CONTROLLER  DESIGN 


Introduction 

This  chapter  details  the  design  of  a  complete  controller  for 
combined  lateral  and  longitudinal  maneuvers  for  the  STOL  like 
flight  condition  (0.15  mach,  sea  level),  along  with  longitudinal 
controllers  for  two  other  flight  conditions  (0.6  mach,  sea  level 
and  0.9  mach,  30,000  ft.).  The  0.15  mach  flight  condition  has  not 
previously  been  studied,  but  the  other  two  have  been  (Ref  2). 
However,  this  thesis  uses  a  different  longitudinal  model  than  any 
previous  efforts  (Ref  2  and  3),  which  led  to  the  redesigning  of 
the  longitudinal  controllers  at  the  0.6  and  0.9  mach  flight 
conditions . 

The  change  in  the  longitudinal  model  results  from  a  forced 
change  at  the  0.15  mach  flight  condition  which  is  then  carried 
over  to  the  other  two  flight  conditions.  This  forced  change 
results  from  a  column  of  zeros  in  the  B  matrix  at  0.15  mach 
corresponding  to  the  maneuver  flaps  (Equation  (48)).  If  this 
column  of  zeros  is  correct,  then  that  means  that  the  maneuver 
flaps  have  no  effect  upon  the  motion  of  the  aircraft  (for  small 
pertubations)  at  this  flight  condition.  Unfortunately,  since  the 
aircraft  does  not  exist,  there  is  no  way  to  verify  this.  Consi- 
dering  that  there  are  at  least  two  other  errors  in  this  same  B 
matrix  (Appendix  A),  this  column  of  zeros  would  have  to  be  consi¬ 
dered  suspect.  This  thesis  is  based  upon  the  idea  of  using  just 
three  of  the  longitudinal  flight  control  surfaces  to  execute  the 
maneuvers,  and  the  decision  to  leave  out  the  maneuver  flaps  was 
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based  upon  this  column  of  zeros.  If  this  zero  column  were  erro¬ 
neous,  the  impact  on  these  results  would  be  minimal. 

Since  the  maneuver  flaps  have  no  apparent  effect  at  0.15 
mach,  they  can  be  eliminated  from  the  longitudinal  model  for  the 
purposes  of  performing  the  maneuvers.  However,  the  maneuver 
flaps  are  needed  to  trim  the  aircraft,  since  without  them  the  low 
angle  of  attack  needed  for  landing  would  not  be  possible.  Based 
on  these  considerations,  it  was  decided  to  model  the  aircraft  as 
using  the  maneuver  flaps  for  trimming  only  and  not  using  them  for 
any  of  the  maneuvers.  Since  this  was  contrary  to  previous 
efforts,  the  0.6  mach  and  0.9  mach  flight  conditions  from 
Bauschlicher's  thesis  were  reaccomplished  using  this  new 
longitudinal  model.  This  would  allow  comparisons  between  the 
different  longitudinal  models,  along  with  drawing  conclusions  as 
to  which  model  generates  the  best  aircraft  performance. 
Unfortunately,  a  design  for  the  2.3  mach  flight  condition  was  not 
found.  Part  of  the  reason  for  this  was  the  limited  time 
available  to  the  author,  but  there  had  to  be  other  reasons  also, 
since  the  other  designs  were  accomplished  in  a  relatively  short 
time  compared  to  the  time  spent  attempting  to  find  a  controller 
for  2.3  mach.  One  of  the  reasons  might  be  the  fact  that  the 
eigenvalues  of  the  system  A  matrix  for  2.3  mach  are  different 
than  the  other  flight  conditions.  As  stated  later  in  this  chap¬ 
ter,  the  short  period  roots  for  all  of  the  flight  conditions 
consist  of  two  real  roots,  one  in  the  left  half  s  plane  and  one 
in  the  right  half  s  plane.  In  fact,  the  short  period  roots  for 
the  2.3  mach  flight  condition  are  very  nearly  identical  to  those 


for  the  0.9  mach  flight  condition.  The  phugoid  roots  for  2.3 
mach,  however,  have  become  real  with  both  lying  in  the  left  half 
s  plane.  This  is  the  only  flight  condition  studied  for  which 
this  occurs. 

The  design  of  a  complete  12  state  combined  lateral  and 
longitudinal  controller  would  be  beyond  the  capabilities  of  MULTI 
(currently  limited  to  10  states),  and  perhaps  the  author  also. 
Consequently,  the  design  is  broken  down  into  two  parts,  the 
longitudinal  and  lateral  modes.  After  the  lateral  and 
longitudinal  controllers  are  designed  and  tested  independently, 
the  lateral  and  longitudinal  models  are  combined  and  a  combined 
controller  is  then  tested  based  on  the  combination  of  the 
separate  controllers.  The  combined  controller  is  tested  using 
one  coupled  maneuver,  to  verify  that  lateral  and  longitudinal 
commands  can  be  executed  simultaneously.  Although  this  phase  is 
simplified  by  the  decoupled  nature  of  the  original  equations  of 
motion  (as  represented  by  the  matrices  in  equations  (48)  and 
(49)),  the  procedure  should  still  be  the  same  if  the  equations 
were  not  decoupled. 

The  computer  program  MULTI  Is  used  extensively  in  the 
design  of  the  controllers  for  this  thesis,  and  is  the  sole  tool 
used  to  test  the  designs  via  simulation.  MULTI  was  originally 
written  by  AFIT  students  in  1981  and  was  based  upon  a  simulation 
program  written  by  Professor  Brain  Porter  (Ref  7).  The  program 
has  undergone  many  alterations  since  then,  including  the  addi¬ 
tion  of  an  option  to  calculate  and  display  the  figures  of  merit, 
and  passing  of  the  input  data  local  filenames  thru  command  line 
arguments  into  the  program.  Both  were  written  by  this  author  and 


detailed  in  Appendix  C.  Appendix  D  lists  the  addition  of  program 
code  that  would  allow  faster  design  times  for  the  experienced 
MULTI  user,  and  would  also  allow  the  true  implementation  of  a  D 
feedforward  matrix.  Although  MULTI  currently  asks  for  a  D  mat¬ 
rix,  the  simulation  does  not  include  a  D  matrix  in  the  calcula¬ 
tions  of  the  output  used  for  feedback.  A  complete  description  of 
MULTI  and  its  options  can  be  obtained  from  the  AFIT  Electrical 
Engineering  Department. 

This  chapter  contains  the  design  and  testing  of  the  separate 
lateral  and  longitudinal  controllers  for  the  0.15  mach,  sea  level 
flight  condition,  along  with  the  longitudinal  controllers  for  the 
0.6  mach,  sea  level  and  0.9  mach,  30,000  ft.  flight  conditions. 
The  approach  used  for  each  flight  condition  is  that  the  control¬ 
ler  designed  will  be  used  for  all  of  the  maneuvers,  instead  of 
trying  to  find  a  different,  possibly  more  optimized  controller 
for  each  maneuver.  Finally,  the  lateral  and  longitudinal  contro¬ 
llers  for  the  0.15  mach  flight  condition  are  combined  and  tested 
with  one  coupled  lateral/longitudinal  maneuver.  Chapter  V  in¬ 
vestigates  using  one  controller  for  all  of  the  flight  condi¬ 
tions  and  the  effect  of  delaying  the  control  input  to  the 
actuators . 

Lateral  Controller 

The  lateral  controller  design  is  accomplished  only  for  the 
first  flight  condition,  0.15  mach,  sea  level.  The  lateral  con¬ 
trollers  for  the  other  flight  conditions  used  in  this  thesis  can 
be  found  in  Bauschlicher ' s  thesis  (Ref  2). 

Manuever  Descriptions.  Four  maneuvers  are  used  in 


Bauschlicher's  thesis  to  test  the  lateral  controller,  the  flat 
turn,  yaw  pointing,  horizontal  translation,  and  the  rollover. 


The  same 

four  maneuvers 

were 

attempted  at  0.15 

mach 

with 

the 

controller 

designed  for 

this 

thesis,  with  the 

result 

that 

the 

rollover  could  not  be  performed  at  all.  This  is  probably  due  to 
the  very  slow  speed  of  the  aircraft  and  subsequent  lack  of 
aerodynamic  forces  exerted  on  the  flight  control  surfaces. 
Consequently,  only  the  first  three  maneuvers,  the  flat  turn,  yaw 
pointing,  and  horizontal  translation,  are  tested  in  this  thesis. 

The  flat  turn  maneuver  is  commanded  by  deciding  on  what  g 
forces  are  desired  in  the  xy  plane  (body  axes),  and  calculating 
the  necessary  yaw  rate  from: 

Ay  -  w**d/2  -  U*r  (50) 

where 

Ay  ■  acceleration  in  xy  plane 
w  ■  r  (body  yaw  rate) 
w*d/2  ■  U  (body  velocity) 

and  with  d/2  used  instead  of  r  for  radius.  The  body  axis  side 
velocity  v  is  commanded  to  zero,  along  with  the  euler  roll  angle 
4>  •  Henceforth,  all  coordinate  systems  are  body  axis,  except  for 
the  euler  angles,  0  (pitch),  Y'(yaw),0(roll). 

Yaw  pointing  is  accomplished  by  deciding  upon  the  angle  of 
pointing  desired,  then  commanding  the  sideslip  angle  p  to  that 
value  via  the  relationship 

P  -  v/U  (51) 

which  is  good  for  small  values  of  sideslip.  Roll  angle  is 
commanded  to  zero,  and  the  command  for  r  is  based  upon  a  curve 


whose  integral  (area  under  curve)  will  be  equal  to  the  sideslip 
angle.  Note  that  an  angle  of  less  than  10  degrees  insures  an 
error  of  less  than  one  percent  in  equation  (51). 

Horizontal  translation  consists  of  commanding  the  side 
velocity  v  depending  upon  desired  g  forces,  while  commanding  both 
roll  angle  and  yaw  rate  to  zero. 

Although  the  rollover  maneuver  is  not  tested  in  this  thesis, 
it  is  performed  by  commanding  both  yaw  rate  and  side  velocity  to 
zero  while  commanding  the  desired  roll  angle. 

Lateral  Model .  These  sections  detailing  the  design  of  the 
lateral  model  and  controller  use  the  flight  condition  of  0.15 
mach  at  sea  level  for  the  numerical  examples  to  illustrate  the 
procedure.  The  same  techniques  are  used  to  develop  the  model  and 
controllers  for  the  other  two  flight  conditions,  with  the  actual 
matrices  used  listed  in  Appendix  D. 

The  first  step  is  to  reduce  the  system  model  from  its 
current  form  of  equation  (48).  The  six  lateral  states  are  p,  r, 
v,  y,  <P ,  and"1}*-.  The  lateral  inputs  are  cf  a ,  cfr,  Svct  F,  ,  , 
CD I,  ,  and  CDI^.  To  reduce  the  number  of  states  without  elimina¬ 
ting  crucial  Information  about  the  aircraft,  the  eigenvalues  can 
be  used.  The  complete  set  of  eigenvalues  for  the  lateral  model 
of  equation  (48)  are: 

-0.1990  +j  0.9563  (dutch  roll  roots) 

-0.6927  (roll  subsidence) 

0.01069  (spiral  divergence) 

+  0.1471  E-06 

The  two  roots  near  zero  contain  no  essential  aircraft  information 
and  so  can  be  eliminated  from  the  model.  Based  on  previous 


experience,  the  choice  is  made  to  eliminate  y  and  from  the 
state  model.  Computing  the  eigenvalues  of  this  new  four  state  A 
matrix  yields: 

-0.1990  +j  0.9563 

-0.6927 

0.01069 

and  the  Important  aircraft  characteristics  have  been  preserved. 

Now  that  the  state  model  has  been  reduced  to  four  states, 
the  inputs  are  going  to  have  to  be  reduced  to  at  most  three  (and 
still  have  the  correct  form  for  an  irregular  design).  This  is 
done  by  observing  that  the  columns  in  the  B  matrix  corresponding 
to  F,  and  F^  have  equal  and  opposite  values  (in  the  lateral  part 
of  the  B  matrix).  This  means  that  as  long  as  the  two  engines 
generate  equal  thrust  these  terms  cancel  out  to  zero,  having  no 
effect  on  the  model.  Likewise,  since  the  terms  corresponding  to 
the  two  coefficients  of  drag  for  the  inlets  are  also  equal  and 
opposite,  they  have  no  effect  with  equal  thrust  in  the  two  en¬ 
gines.  This  reduces  the  Inputs  to  three  for  the  lateral  model, 
&a,  dr,  and  6v c.  This  model  is  satisfactory  as  long  as  the  two 
engines  produce  the  same  thrust.  Imposing  this  constraint  would 
require  extra  circuitry  in  addition  to  the  controller.  Quite 
possibly,  this  constraint  on  the  thrust  would  be  valid  even 
without  any  additional  engine  control. 

The  choice  of  the  three  outputs,  v,  c fi,  and  r,  is  dictated 
by  Bauschlicher's  thesis,  since  one  of  the  purposes  of  this 
thesis  is  to  compare  results  with  Bauschlicher .  This  choice  of 
outputs  satisfies  the  requirement  for  the  number  of  inputs  to 


equal  the  number  of  outputs.  These  outputs  result  in  completely 
different  output  matrices  than  those  in  equation  (49).  Since 
there  are  no  accelerations  in  the  output,  there  is  no  feedforward 
matrix,  and  the  C  matrix  is  also  different  from  equation  (49). 

These  changes  result  in  a  new  set  of  matrices  describing  the 


linear  dynamics  of  the  aircraft  and  are  given  by: 


31.6  -.172  -164  31.5 


.00415  -.212  .0654 


-.00818  .364  -.696 


.111  .0363 


.00114  -.0112  .0022 


.0122  .0118  -.00212 


(52) 


0  10  0  <f> 


10  0  0 


0  0  10 


|pj  (53) 
The  above  set  of  equations  represents  the  lateral  model  for  the 
FPCC  aircraft  at  0.15  mach  at  sea  level.  This  system  is  checked 
for  decoupling  zeros,  and  since  there  aren't  any  it  is  both 
controllable  and  observable,  two  necessary  requirements. 
Although  Bauschlicher  stated  that  the  transmission  zeros  must 
also  lie  in  the  left  half  plane,  this  test  is  performed  on  the 
system  with  F  ("[C(  +MA((  HC^+MA^  ])  substituted  for  C.  Making 
this  substitution  and  then  checking  for  transmission  zeros  yields 
none  in  the  right  half  plane,  another  requirement.  (Although  it 


may  be  possible  to  design  with  transmission  zeros  at  the  origin) 
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Since  the  first  Markov  parameter  (see  Chapter  II)  CB  ■  C  B 
does  not  have  full  rank  of  1,  the  design  is  irregular.  Conseque¬ 
ntly,  an  M  matrix  must  be  found  which  meets  the  requirements 
specified  in  Chapter  II.  This  measurement  matrix  supplies 
information  to  the  feedback  loop  on  the  derivative  of  <p ,  which  is 
not  necessarily  equal  to  the  roll  rate.  The  following  section 
details  the  development  of  the  M  matrix. 


Measurement 

Matrix  Determination. 

The 

first  step 

i  s 

partitioning  the 

matrices  of  equations 

(52)  and 

(53)  so 

the 

methods  of  Chapter 

II  can  be  used  to  find 

the  B* 

matrix. 

The 

partitioned  matrices  used  in  equation  (43)  are: 


C  -  Cf 


—  m  —a 


B  -  A 
—  ~ia 


0 
1 
0 

•  — 

Coj 

[0  0.19  1] 

The  resulting  B*  matrix  is  found  to  be: 
10  0  0i 

B* 


(54) 

(55) 

(56) 


(57) 


0  0.19  1 

0  0  0 

with  { d,  ,  dA ,  d^  , }  »  (3,0,3}  as  used  in  equation  (43). 

Next,  the  measurement  matrix  M  must  be  found,  so  that  can 
be  calculated.  Since  the  dimensions  of  M  are  lx(n-l),  which  for 
this  system  is  3x1,  M  can  be  generalized  as: 

m, 

M  -  |m,  (58) 
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Using  F 

-a. 


5a 


yields 


-Fa 


C  +MA  ,  and 

—  i  —  i  j. 


0  0 
0  0 
1  0 


0. 19m, 
0 . l9mA 


(59) 


m . 


1+0.1 9m^  m3 


(60) 


Since  the  rank  of  F  must  be  equal  to  lt  or  three,  and  the 
allowable  non  zero  elements  of  F^  (due  to  M)  must  have  a 
corresponding  non  zero  element  in  B*,  the  choice  of  M  must  be: 

0 

M  ■  I m,  (61) 


which  yields 
1 


F 

"A. 


0 

0.  19m 


a  “a 


(62) 


0 

0  1  0 

Using  the  above  information  the  asymptotic  transfer  function  from 
equation  (38)  can  be  derived: 


-i 

F 

-A 


~1  0  0 
0  0  1 
0  1/m  -0.19 


-i 

C  F 
-a- a. 


fl  0  0 
0  0  0 
0  0  1 


which  is  in  the  desired  diagonal  form. 


F(  -  C(  +  MA  -  C 


Since  A 

—it 


C  -C  F~'  F 

“I  — A-a.  —i 


Cta  a  f"'  3  -Ct  (0  l/ma  o)] 

£'  1,  1  -U-l+T/m,af ' 

•[l/CA-i+T/m^)  ] 


which  yields: 


[C  -C  f'  F][>I  -I  -TA  +TA  f"  F  3  [TA  F~'l 
-•  -•»-*  -|J<-  -A-i  -/(  -ia  -a  -/J  L  -u  -a  J 


0  0 


0  T/ma  0  Cl/(X-l+T/ma)3 


0  0 


This  has  the  desired  diagonal  form. 


<r,  0  0 


F  B  K  -  0  <r,  0 

-A-A-o  J 


0  0  <r, 


[U  -I  +F  B  K  1 

-  X.  -JJ  "A  -  J  -  0J 


A  -1+07  0 


A  -1  +  a^  0 


A  -1+cr 


i/C\-i+cr;)  0 


i/(A-i  +  <r,)  0 


(63) 


1/0-1+05 ) 


<r,  /(A-i+o: )  0  0 


[C  F  1  [A  I  -I  +F  B  K  3  [F  B  K  3-  0 


0  0 


0  <£/u-i+cp 
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which  is  also  in  the  desired  diagonal  form. 

Now,  since  ^  *  P  +  T 

["^/(A-I+gt)  0  0 

0  -  0  (T/oa)/(A  -l+T/mA)  0 

0  0  <Tj  /(A  -1  +  0} ) 


(64) 


The  asymptotic  transfer  function  in  equation  (64)  indicates 
that  the  desired  design  goal  of  output  decoupling  with  infinite 
gain  has  been  achieved.  Note  that  a  decrease  in  sampling  time 
corresponds  to  an  increase  in  gain,  so  that  saying  the  sampling 
time  goes  to  zero  is  equivalent  to  stating  that  the  gain  goes  to 
infinity.  From  equation  (64)  the  closed  loop  eigenvalues  for 
this  system  are: 

A,-  i-g; 

Aa»  1-T/m* 

l-13} 

These  eigenvalues  indicate  that  the  time  response  of  the  outputs 
are  greatly  affected  by  the  values  of  <T,  ,  m2,  and  . 

Going  back  to  an  intermediate  step  in  the  derivation  of 
equation  (64),  the  one  (n-1)  transmission  zero  for  this  system  is 
(equation  (63)): 

Zt  -  1-T/m^ 

Controller  Design.  Detailed  below  is  the  basic  algorithm 
used  for  the  controller  design  of  all  of  the  controllers  designed 
for  this  thesis.  This  procedure  is  based  mostly  on  the  control¬ 
lers  within  this  thesis  and  as  such  should  not  be  considered  well 
tested.  However,  there  is  no  reason  not  to  expect  this  technique 


to  help  the  designer  get  started  with  a  design.  Implementing 
this  process  depends  upon  available  design  tools,  and  the  one 
used  for  this  thesis  is  the  interactive  program  MULTI.  Without 
MULTI  this  thesis  would  not  have  been  possible.  The  design 
process  follows  these  steps: 

1.  Set  all  of  the  design  parameters  equal  to  one:  M  ele¬ 
ments;^.  elements  ;  «■ ;  and  £  . 

2.  Calculate  Kq  and  K  .  Reduce  until  the  largest  order 
of  magnitude  of  elements  in  either  controller  matrix  is  10E+03. 

3.  Run  a  simulation  and  use  option  #28  in  MULTI  (figures  of 
merit)  to  check  for  a  bounded  (stable)  response.  Continue 
reducing  £  ,but  by  relatively  small  Increments,  until  the 
response  is  stable.  Experience  will  best  answer  the  question  of 
what  relatively  small  increments  are,  but  a  good  starting 
guideline  would  be  factors  of  two.  (Based  upon  the  designs  in 
this  thesis,  step  three  will  reduce  the  largest  order  of  magni¬ 
tude  to  10E-H)2,  but  information  obtained  from  other  thesis  ef¬ 
forts  suggests  the  above  guideline  of  10E+03  in  step  two) 

4.  Further  mold  the  controller  matrices  using  the  elements 
of  2:  .  Increase  or  decrease  the  individual  elements  so  as  to 
Increase  or  decrease  the  corresponding  columns  in  the  controller 
matrices  by  the  same  factor  until  the  orders  of  magnitude  lie 
within  10E-02  to  10E+02  (or  as  close  as  possible).  Generally, 
numbers  that  are  smaller  than  the  suggested  values  are  more 
desirable  than  numbers  that  are  larger  than  suggested. 

5.  Run  a  simulation  to  observe  which  outputs  need  improved 
response.  Generally,  increasing  an  element  in  the  !|  matrix  will 
increase  the  speed  of  response  of  some  of  the  outputs,  but  with 


the  resulting  problem  of  increased  overshoot.  Experimentation 


will  show  which  S  element  affects  which  output.  Repeat  this 
process  until  the  best  time  responses  are  achieved. 

6.  Experiment  with  the  value  of  while  observing  output 
response.  In  this  thesis  increasing  sometimes  reduced  over¬ 
shoot  on  the  control  surface  deflections  with  little  degredation 
in  output  response,  but  sometimes  this  parameter  had  little 
notlcable  effect. 

7.  Now  vary  the  M  elements  to  see  if  they  will  improve 
output  response.  Generally  the  M  elements  directly  affect  the 
corresponding  output  response  depending  upon  which  column  they're 
in.  For  example,  decreasing  an  entry  in  the  second  column  might 
Increase  the  speed  of  response  of  the  state  whose  derivative  is 
measured  by  that  column.  Host  designs,  including  those  from 
other  theses,  have  values  for  M  elements  within  the  range  of  0.25 
to  1.0. 

8.  Repeat  steps  5  thru  7  until  the  best  responses  are 
achieved,  keeping  in  mind  that  tradeoffs  will  always  be  made 
between  desirable  and  undesirable  output  responses  and  control 
surface  deflections.  Even  though  option  #27  in  MULTI  does  limit 
the  control  surface  deflections  based  on  physical  limitations,  a 
design  that  reduces  the  number  of  times  a  surface  "hits  the 
stops”  1 8  considered  superior  to  one  that  doesn't. 

Obviously  different  plants  are  going  to  respond  in  different 
ways  to  changes  in  the  design  parameters.  This  means  that  a 
satisfactory  design  is  more  the  result  of  trial  and  error  than 
any  "cookbook"  approach.  One  conclusion  tentatively  drawn  from 


current  and  past  theses  is  that  fighter  aircraft  seem  to  respond 
to  design  parameter  changes  similarily,  as  do  transport  aircraft, 
although  fighters  differ  from  transports. 

Since  the  asymptotic  transfer  function  is  valid  only  for 
values  of  gain  approaching  infinity,  for  values  of  gain  much  less 
than  infinity  there  will  be  coupling  between  the  outputs.  This 
is  probably  one  of  the  main  causes  of  frustration  that  results 
when  trying  to  shape  one  output  response  only  to  find  that 
another  is  severely  degraded.  On  real  aircraft,  the  value  for 
the  sampling  time  will  be  dictated  by  the  performance  of  the 
processor  (computer)  used.  Unfortunately,  it  may  be  impossible 
to  run  the  processor  at  the  speed  that  yielded  those  good  results 
when  designing  with  MULTI. 

For  the  lateral  controller  different  controller  matrices 
yielded  the  best  results  for  the  different  commanded  maneuvers  at 
this  one  flight  condition.  Since  the  philosophy  of  this  thesis 
is  to  try  and  follow  Bauschl icher ' s  approach  so  as  to  make 
comparisons,  one  set  of  controller  matrices  Is  chosen  for  each 
flight  condition.  This  results  in  compromise  among  the 
individual  responses  to  each  commanded  maneuver.  The  choice  of 
which  controller  matrices  to  use  is  subjective,  based  on  the 
author's  evaluation  of  what  constitutes  the  best  "overall" 
response . 

The  final  values  selected  for  the  lateral  controller  for  the 
0.15  mach  at  sea  level  flight  condition  are  listed  below. 
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Figures  5  through  11  show  the  time  responses  for  the  outputs 
and  the  flight  control  surface  deflections  for  the  0.15  mach,  sea 
level  flight  condition.  The  figures  of  merit  and  commanded 
inputs  for  each  maneuver  are  listed  in  Table  VII. 
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TABLE  VII 

SUMMARY  OF  OUTPUT  RESPONSES  FOR  LATERAL  CONTROLLER 


Maneuver  Command 


Overshoot 


Rise  Time 


0.15  Mach 

Yaw 

1,-11.5,20,20 

none 

1  sec 

0  ft. 

Pointing 

0,0, 0,0 

-0.001  rad 

.2, .0524, 1,  1.2 

8.2% 

0.27  sec 

Flat 

0,0, 0,0 

-38.9  ft/sec 

T.rn 

0,0, 0,0 

.068  rad 

1, .195,20,20 

none 

* 

Ho riz  . 

1,48.3,20,20 

none 

10  sec 

Transl  . 

0,0, 0,0 

.052  rad 

0,0, 0,0 

-.12  rad/sec 

Note:  The  command  vector  Is  [  v  <p  r  ]. 

Overshoot  is  the  first  peak  for  non  zero  commanded  outputs 
or  is  the  greatest  deviation  from  zero  for  zero  commanded 
outputs.  Rise  time  is  from  initial  value  to  90%  of 
commanded  value. 

*  For  this  maneuver  the  output  did  not  reach  90%  of  the  command. 


3  DEGREE  YAH  POINTING  (0.  0.15  NACH.  0  FT) 

Figure  11.  3  Degree  Yaw  Pointing  (0.15  Mach) 
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Longitudinal  Controller 

Maneuver  Description.  Continuing  with  the  concept  of  direct 
comparison  with  Bauschlicher's  results,  the  same  three 
longitudinal  maneuvers  that  he  used  are  chosen  for  this  thesis. 
These  maneuvers  are  the  pitch  pointing,  vertical  translation,  and 
direct  lift.  These  maneuvers  are  typically  claimed  to  be  very 

I 

useful  in  the  combat  arena,  both  for  air-to-air  and  air-to-ground 
combat . 

Pitch  pointing  involves  equal  changes  in  the  angle  of  attack 
and  the  pitch  angle,  resulting  in  the  flight  path  angle  -oc) 
remaining  constant.  Therefore,  the  nose  pitches  while  the 
aircraft  remains  on  the  current  trajectory.  Note  that  using  the 
euler  pitch  angle  implies  that  the  aircraft  has  zero  initial  roll 
angle.  To  generalize  this  maneuver  for  any  initial  attitude,  the 
integral  of  pitch  rate  could  be  used  instead  of d  .  The  change  in 
velocity  1 8  commanded  to  zero. 

Vertical  translation  results  in  the  aircraft  moving  in  the 
vertical  direction  (along  z  axis)  with  no  change  in  the  pitch 
angle.  Consequently,  the  pitch  angle  change  is  commanded  to 
zero,  as  is  the  change  in  velocity.  The  angle  of  attack  is 
commanded  so  as  to  actually  command  a  change  in  w,  the  velocity 
in  the  ?  direction.  Since 

<*  -  w/U  (71) 
for  small  angles  of  attack,  changes  in  °<  equate  to  changes  in  w 
(and  thus  commanded  g  forces  in  the  z  direction). 

Previous  thesis  students  have  called  the  third  maneuver  a 
direct  lift,  even  though  direct  lift  has  typically  been  used  to 


describe  a  different  maneuver  while  the  constant  g  pull  up  has 


been  used  to  describe  this  last  maneuver.  This  thesis  will 
follow  previous  theses  in  calling  this  maneuver  the  direct  lift. 
To  command  this  maneuver,  equation  (50)  can  be  rewritten  as 

Az  -  q*U  (72) 
from  which  the  pitch  angle  command  necessary  to  generate  the 
desired  acceleration  in  the  z  direction  (Az)  can  be  found.  Both 
the  change  in  u  and  angle  of  attack  are  commanded  to  zero,  so 
that  the  flight  path  angle  follows  the  pitch  angle. 

Longitudinal  Model .  Just  as  with  the  lateral  controller, 
modifications  to  the  data  generated  by  the  FPCC  simulation 
program  are  required.  The  modifications  made  to  equations  (48) 
and  (49)  in  this  chapter  apply  to  the  0.15  mach  flight  condition, 
but  the  same  changes  are  also  made  to  the  data  for  the  other 
flight  conditions  (Appendix  D). 

The  longitudinal  states  used  in  equation  (48)  are  q,  w,  u, 
0,  h,  and  x.  Just  as  with  the  lateral  model,  an  analysis  of  the 
eigenvalues  of  the  A  matrix  can  be  used  to  reduce  the  system 
state  model.  The  eigenvalues  of  the  original  six  state  A  matrix 
are : 

-0.02010  +  j  0.2034  (phugoid  roots) 

1.166  (unstable  short  period  root) 

-1.887  (stable  short  period  root) 

0.0001981 

-0.001131 

Since  this  aircraft  is  longitudinally  unstable,  one  of  the  short 
period  roots  has  migrated  over  into  the  right  half  s  plane.  The 
two  roots  near  zero  supply  no  essential  information  and  so  can  be 
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removed  from  the  model. 


Once  again  from  experience  the  decision 


is  made  to  eliminate  the  two  states  h  and  x  yielding  the  follo¬ 
wing  eigenvalues: 

-0.02059  +  j  0.2040 
1.166 
-1.887 

So  the  plant  model  has  been  reduced  to  four  state,  q,  w,  u,  and 
9  with  no  loss  of  essential  aircraft  characteristics. 

From  equation  (49)  there  are  seven  longitudinal  control 
inputs,  Sc.,  cfj ,  Saf ,  F,  ,  FA  ,  CDI,  ,  and  CDI^.  However,  for  this 
particular  flight  condition  there  is  a  column  of  zeros  in  the  B 
matrix  corresponding  to  the  maneuver  flaps,  implying  they  have  no 
effect  upon  the  aircraft  in  this  configuration.  Of  the  remaining 
six  control  inputs,  the  two  engine  thrusts  can  be  combined  into 

9 

one  input  called  Fror  ,  by  simply  adding  those  two  columns  of  the 
B  matrix  together.  Once  again  the  contributions  of  the  terms  in 
the  equations  due  to  CDI,  and  CDI^  are  minimal  and  can  be 
dropped.  This  leaves  just  three  control  inputs,  Sc,  Sj,  and  FTJr> 
Since  there  are  three  Inputs,  there  can  be  at  most  three 
outputs.  The  outputs  chosen  are  the  same  as  Bauschlicher  used 
except  for  thrust.  These  outputs  are  9,  u,  and  m.  Since  angle 
of  attack  isn't  one  of  the  states,  it  is  generated  as  an  output 
by  assuming  that  <x  is  equal  to  w/U.  This  means  that  in  each  C 
matrix  used  for  this  thesis  this  entry  is  the  reciprocal  of  the 
airspeed  at  trim. 

The  relationship  between  the  input  command  to  the  thrust  and 
sV  the  engine  thrust  is  modeled  as  an  actuator  in  this  thesis.  This 
is  preferable  to  imbedding  the  engine  thrust  into  the  plant  A 
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matrix.  Imbedding  the  thrust  (or  any  other  actuator)  means  that 
the  model  assumes  instantaneous  aircraft  response  to  the  change 
in  thrust  (or  flight  control  surface),  whereas  keeping  the 
actuator  dynamics  separate  and  using  option  if 4  in  MULTI  means  the 
model  incorporates  a  (more  realistic)  one  sampling  period  time 
delay  between  actuator  movement  and  aircraft  response. 

With  these  new  states,  inputs,  and  outputs  the  longitudinal 
model  used  in  this  thesis  is  expressed  by  the  following 
equa  tions : 
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0 

0 
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-31.6 

-.0214 

.0867 

-31.5 

-6.01 

-.215 
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(74) 


An  interesting  aside  is  that  the  system  above  yielded  no  transmi¬ 
ssion  zeros  and  no  decoupling  zeros. 

The  above  equations  represent  the  longitudinal  model  for 
this  aircraft  at  0.15  mach,  sea  level.  This  system  is  checked 


(with  F  substituted  for  C)  for  transmission  zeros  and  none  are 
found  in  the  right  half  plane.  No  decoupling  zeros  are  found, 


Indicating  that  the  system  is  controllable  and  observable. 


Since  the  matrix  C  B  does  not  have  full  rank  1,  the  system 
is  irregular  and  the  measurement  matrix  must  be  found.  The 
concept  behind  the  ir~egular  design  is  explained  in  Chapter  II. 

Measurement  Matrix  Development .  Since  n-1  is  one,  the  M 
matrix  is  once  again  3x1,  this  time  measuring  the  derivative  of 
the  pitch  angle.  The  first  step  is  to  find  B* ,  which  can  be 
derived  from: 


C 
-  r 


A 

—  tt 


1 
0 
0 

-  Co] 


(75) 


A  -[001] 
-ix 


(76) 

(77) 


yielding 


B* 


with 


0 

0 
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0 

0 
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1 

0 
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{ d,  ,  da,  d3  }  -  {0,  3,  3} 
Using  the  general  form  for  M: 
m, 


(79) 


M 


and 
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00607  m. 


For  F  to  have  full  rank,  yet  be  as  sparse  as  possible,  and  using 
the  guidelines  of  Chapter  II  in  choosing  the  M  elements  based 
upon  B* ,  the  obvious  choice  for  M  is: 


(83) 


which  yields 


0  0 


1  0 


(84) 


0  .00607  0 


1  0 


0  1/. 00607 


1/m,  0  0 


(85) 


The  asymptotic  transfer  function  matrix  can  be  derived  as  in 
the  lateral  model  earlier  in  this  chapter  and  is  found  to  be: 


(T/m,  )/(X-l+T/m,  )  0 


VOv-l  +  o^)  0 
0  <5/0-1+05) 


(86) 


From  equation  (86)  the  transmission  zero  is 
Z*  ■  1-T/m, 


(87) 


Equation  (86)  shows  that  the  outputs  become  decoupled  as  the  gain 
goes  to  infinity  (or  sampling  time  T  goes  to  zero).  Also,  m,,o~, 
and  have  the  greatest  Impact  on  the  output  time  responses. 
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Controller  Design .  The  same  algorithm  listed  under  lateral 
controller  design  is  used  for  the  longitudinal  controller  design. 
Actually,  this  algorithm  is  used  for  every  controller  design  for 
this  thesis.  The  choice  of  one  set  of  controller  matrices  for 
all  three  maneuvers  for  each  flight  condition  requires  compro¬ 
mises  in  the  performance  of  each  maneuver.  The  design  values 
used  are  listed  below. 
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0.1,  T 
0.25 
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0 


0.01 
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16.856 

.00589 

-18.518 

-.  564 

427.964 

89.785 

“8.428 

.00295 

-9.259 

-.282 

213.982 

44.893 

-1.340 

-5.0981 

-100.0318 

-.670 

-2.549 

-50.0159 


(88) 

(89) 

(90) 

(91) 


(92) 


(93) 


Appendix  D  lists  all  of  the  design  parameters  for  the  other 
flight  conditions. 

Figures  12  through  20  show  the  time  responses  for  the 
outputs  and  control  inputs  for  the  three  longitudinal  maneuvers 
at  0.15  mach,  sea  level.  Figures  21  through  39  show  the  same 


time  responses  for  the  same  maneuvers  at  the  other  two  flight 


1.75  DEGREE  PITCH  POINTING  (0.15  MACH.  0  FT) 


1.75  DEGREE  PITCH  POINTING  (0.15  MRCH .  G  FT) 

Figure  12.  1.75  Degree  Pitch  Pointing  (0.15  Mach) 
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0.5  0  VERTICAL  TRANSLATION  10.15  MACH.  0  FT) 


Figure  18.  0.5  g  Vercical  Translation  (0.15  Mach) 
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conditions,  0.6  mach,  sea  level,  and  0.9  mach,  30,000  ft.  Table 
VIII  shows  the  figures  of  merit  from  these  time  responses  along 
with  the  commanded  inputs  for  each  maneuver. 


Controller  Response  Analysis 

Lateral  Controller  Response  Analysis .  Figure  5  shows  the 
side  velocity  response  for  the  horizontal  translation  maneuver, 
the  only  non-zero  commanded  input.  Obviously  the  response  does 
not  follow  the  command,  either  in  the  transient  stage  or  in  the 
steady  state  (defined  here  as  ten  seconds  after  initial  command 
input).  Experimenting  with  other  commanded  g  levels  resulted  in 
similar  ti»e  responses  of  the  side  velocity.  Optimizing  the 
controller  for  just  this  one  maneuver  did  not  significantly 
Improve  the  response.  This  may  be  due  to  the  presence  of  a  slow 
mode,  due  to  a  transmission  zero,  in  the  output. 

Figure  7  shows  the  time  response  of  the  yaw  rate  for  the 
flat  turn,  the  only  non-zero  commanded  input.  For  this  flight 
condition,  this  is  the  worst  responding  maneuver  of  the  three. 
Apparently  the  aircraft  tracked  the  commanded  yaw  rate  fairly 
well  until  one  second,  when  the  command  was  leveled  off  to  a 
constant  value.  At  this  time  the  yaw  rate  falls  off  rapidly  while 
the  roll  angle  increases  as  does  the  side  velocity,  until  in  the 
steady  state  the  aircraft  is  turning  and  sliding  sideways  through 
the  alrstream. 


Figure  8  shows  the  control  surface  deflections,  indicating 
that  the  control  surfaces  are  being  taxed  due  to  the  slow 


airspeed 


The  vertical  canard  runs  immediately  to  full 


deflection.  and  the  ailerons  are  run  from  full  deflection  in  one 


direction  to  full  deflection  in  the  other  in  just  one  second. 


Longitudinal  Controller  Response  Analysis .  For  the  0.15  mach 
flight  condition,  Figure  12  shows  the  output  responses  for  the 
1.75  degree  pitch  pointing  maneuver.  Although  the  transient 
response  is  fairly  acceptable,  the  steady  state  response  is 
probably  not,  with  alpha  less  than  theta  in  magnitude.  This 
results  in  a  maximum  flight  path  angle  of  about  0.9  degrees 
Instead  of  the  commanded  zero  degrees. 

Figure  13  shows  the  oscillations  in  thrust  and  jet  flap 
deflection,  which  have  the  same  frequency  as  the  oscillations  in 
forward  velocity.  Here  is  the  first  indication  that  the  jet  flap 
1 8  not  going  to  work  as  well  as  the  maneuver  flap  did  for 
Bauschlicher  in  performing  as  an  elevator.  The  problem  is  that 
deflecting  the  jet  flap  to  Induce  a  pitching  moment  also  deflects 
the  thrust  and  affects  the  forward  velocity.  The  horizontal 
canard  must  deflect  both  to  control  any  pitching  moment  from  jet 
flap  deflection  and  also  to  control  changes  in  lift  due  to 
deflection  of  thrust  (and  consequent  changes  in  velocity).  The 
deflection  rates  of  the  horizontal  canards  and  jet  flaps  present 
the  least  desirable  (and  possibly  obtainable)  aspects  of  this 
maneuver . 

Figure  16  indicates  that  the  aircraft  response  to  the  direct 
lift  maneuver  followed  the  commanded  inputs  quite  well.  However, 
the  0.5  g  commanded  was  about  the  maximum  possible  using  this 
controller.  Even  optimizing  the  controller  for  this  maneuver  did 
not  add  significantly  to  the  aircraft  performance.  The  apparent 
instability  at  the  end  of  the  simulation  time  probably  results 
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from  either  wind  up  or  from  an  induced  instability  in  the 
controller  due  to  a  lag  in  the  integration  routine.  Appendix  C 
contains  more  detail  concerning  this  problem. 

Another  possible  problem  is  the  initial  jet  flap  deflection 
needed  to  trim  the  aircraft  at  0.15  mach,  sea  level.  Obviously 
this  adds  complications  not  present  with  zero  initial  jet  flap 
deflection,  such  as  in  Bauschlicher' s  thesis. 

Figure  18  shows  the  output  responses  for  the  v  tlcal  tran¬ 
slation  maneuver  at  0.15  mach,  sea  level.  The  out  ~s  result  in 
a  lower  g  rating  than  commanded,  but  it  is  stil  a  vertical 


translation  with  almost  no  change  in  pitch  angle 


ace  again  an 


instability  occurs  at  the  end  of  the  simulation.  Appendix  C 
contains  more  Information  on  the  causes  and  cures  for  this 
instability. 

From  the  pitch  pointing  and  vertical  translation  maneuvers, 
one  conclusion  is  that  the  alpha  response  does  not  track  commands 
as  well  as  does  the  theta  response. 


Combined  Lateral  and  Longitudinal  Controller 


Development .  The  last  section  of  this  chapter  shows  how  the 
longitudinal  and  lateral  models  can  be  combined  along  with  the 
lateral  and  longitudinal  controllers.  Since  the  system  A  matrix 
for  the  0.15  mach,  sea  level  flight  condition  is  decoupled 
(equation  (48)),  as  is  the  B  matrix  after  dropping  the  columns 
corresponding  to  CD1,  and  CDIa,  combining  the  separate  models 
together  poses  no  problems.  The  original  columns  for  each  thrust 
input  in  the  B  matrix  are  added  together,  just  as  for  the 
longitudinal  model.  Combining  the  lateral  and  longitudinal 
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system  equations  into  one  set  of  system  equations  yields: 
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The  combined  system  shown  in  equation  (94)  is 
observable  and  controllable,  and  has  no  transmission  zeros  in  the 
right  half  s  plane.  After  substituting  the  F  matrix  for  the  C 
matrix  the  transmission  zeros  are  found  to  be  in  the  left  half 
plane  and  are: 

Z4  -  {1-T/m,  ,1-T/m^  }  (95) 

with 


M 


.25  0 

0  0 

0  0 

0  0 

0  .25 


0 


(96) 


The  r,  2  ,  Kq,  and  K(  matrices  are  just  the  block  diagonal 
combinations  of  the  corresponding  matrices  for  the  individual 
designs  with  the  off  block  diagonal  terms  equal  to  zero. 
However,  the  scaler  has  now  become  a  matrix  and  is  given  by: 


(97) 


To  demonstrate  the  decoupling  of  the  combined  system,  the 
separate  lateral  Inputs  for  the  horizontal  translation  and  the 
longitudinal  inputs  for  the  vertical  translation  are  combined  to 
command  a  jink  maneuver.  As  can  be  seen  in  Figures  40  through 
43,  the  responses  to  the  combined  jink  maneuver  are  nothing  more 
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ROL  SURFACE  INPUTS 

-60.00  -40.00  -20.01 


TABLE  VIII 


SUMMARY  OF  OUTPUT  RESPONSES  FOR  LONGITUDINAL  CONTROLLERS 


Fit.  Cond. 

Maneuver 

Command,  V 

Overshoot 

Rise  Time 

0.15  Mach 

Pitch 

.5, .0305,20,20 

4Z 

.5  sec 

0  ft. 

Pointing 

0,0, 0,0 

-.385  ft/sec 

.5, .0305,20,20 

3Z 

.75  sec 

Direct 

10, .976,20,20 

none 

9.9  sec 

Lift 

0,0, 0,0 

-.861  ft/sec 

0,0, 0,0 

.07  rad 

Vertical 

0,0, 0,0 

.00 9  rad 

Tr ansi . 

0,0, 0,0 

.323  ft/sec 

1, .0977,20,20 

* 

9  sec 

0.6  Mach 

Pitch 

.5, .0305,20,20 

48Z 

.8  sec 

0  ft. 

Pointing 

0,0, 0,0 

.348  ft/sec 

.5, .0305,20,20 

none 

** 

Direct 

10, .865,20,20 

none 

** 

Lift 

0,0, 0,0 

-.324 

0,0, 0,0 

.0132  rad 

Vertical 

0,0, 0,0 

.0239  rad 

Transi . 

0,0, 0,0 

.339  ft/sec 

1, .0384,20,20 

none 

9.1  sec 

0.9  Mach 

Pitch 

.5, .0524,20,20 

12. 2Z 

1.5  sec 

30,000  ft. 

Pointing 

0,0, 0,0 

-1.39  ft/sec 

.5, .0524,20,20 

none 

** 

Direct 

10, .865,20,20 

none 

10  sec 

Lift 

0,0, 0,0 

-3.74 

0,0, 0,0 

.0416  rad 

Vertical 

0,0, 0,0 

-.00459  rad 

Transi . 

0,0, 0,0 

-.283  ft/sec 

1,-. 036, 20, 20 

none 

** 

Note:  Overshoot  is  the  first  peak  for  non  zero  commanded  outputs 

and  1 8  the  greatest  deviation  from  zero  for  zero  commanded 
outputs.  Rise  time  is  from  initial  value  to  90Z  of  the 
commanded  value.  T 

Command  vector  is  [9  u  o<  ]  . 

*  For  this  maneuver  the  shape  of  the  curve  did  not  yield  an 
obvious  value  for  overshoot. 

**  For  this  maneuver  the  output  did  not  reach  90Z  of  the 
commanded  value. 


in 


•.  v.\  s  .*• 


than  the  responses  to  the  vertical  and  horizontal  translation 


commands  (see  Figures  5  through  11  and  12  through  20).  However, 
if  the  aircraft  were  not  decoupled,  then  this  would  not  be  true 
and  the  response  to  a  coupled  maneuver  would  be  different  than 
that  obtained  by  commanding  the  lateral  and  longitudinal  Inputs 
separately.  In  this  case  the  controller  may  have  to  be  modified 
to  compensate  for  coupling  in  the  system  matrices. 


Conclusions 

Not  unexpectedly,  the  aircraft  response  to  the  commanded 
inputs  is  generally  better  for  the  lateral  design  than  for  the 
longitudinal  designs.  This  may  be  due  to  the  use  of  the  jet 
flaps  to  induce  a  pitching  moment  (instead  of  maneuver  flaps) 
which  also  vectors  the  thrust.  Obviously,  having  one  control 
surface  affect  two  control  inputs  is  not  as  efficient  as 
separating  the  control  inputs.  Comparing  these  results  to 
Bauschllcher's  (using  maneuver  flaps),  the  conclusion  is  that 
using  the  maneuver  flaps  for  longitudinal  maneuvers  yields 
significantly  Increased  aircraft  performance. 

Designing  one  controller  laterally  and  longitudinally  for 
each  flight  condition  that  functioned  nearly  as  well  as  separate 
controllers  designed  for  each  maneuver  further  demonstrates  the 
power  of  the  Porter  method,  and  is  an  improvement  over  the 
current  practice  of  gain  scheduling  by  maneuver.  Extending  this 
concept  further  to  one  controller  for  multiple  flight  conditions 
is  examined  in  the  following  chapter. 


CHAPTER  V 


ROBUSTNESS  TESTING  OF  LONGITUDINAL  CONTROLLER 


COMPENSATION  FOR  EFFECTS  OF  CONTROLLER  DELAY 


Introduction 


This  chapter  details  the  design  end  testing  of  a  "universal 
controller"  that  can  be  used  for  nultlple  flight  conditions.  The 
obvious  advantage  to  such  a  robust  controller  Is  that  the 
controller  natrlces  gains  do  not  have  to  be  changed  as  often  as 
If  this  robustness  was  not  possible,  a  practice  called  gain 
scheduling.  Any  valid  claiaa  for  robustness  using  the  Porter 
method  would  have  to  be  substantiated  with  more  complete 
robustness  testing  than  is  possible  for  this  thesis.  The  attempt 
here  is  just  to  demonstrate  that  some  robustness  does  exist  and 
how  to  test  robustness.  This  chapter  also  details  the  effects  of 
a  one  sampling  period  time  delay  In  the  output  of  the  PI 
controller.  This  time  delay  represents  the  time  delay  that  would 
actually  exist  In  the  hardware  realisation  of  the  design  due  to 
the  sample  and  hold  devices  (A/D  converters)  and  the  processing 
delay.  Another  aspect  of  gain  reduction  is  finding  certain  gains 
within  the  controller  matrices  that  can  be  set  to  zero  for 
multiple  flight  conditions.  Unfortunately,  this  approach  did  not 
yield  the  hoped  for  results  with  this  aircraft,  but  the  limited 
result*,  achieved  are  summarized  in  this  chapter. 


Universal  Controller 


The  approach  used  for  this  thesis  in  searching  for  a  universal 
controller  was  to  test  the  controllers  for  each  flight  condition 


this  did  not 


at  the  other  flight  conditions.  Unfortunately, 
result  in  any  one  controller  that  performed  satisfactorily  for 
all  three  flight  conditions.  (Since  only  multiple  longitudinal 
controllers  were  designed,  this  chapter  involves  only  longitudi¬ 
nal  controllers)  Further  testing  with  redesigned  controllers 
also  did  not  result  In  any  one  controller  showing  acceptable 
results  for  all  three  flight  conditions.  Consequently,  the  re¬ 
quirement  for  a  universal  controller  was  downgraded  to  having 
good  responses  for  two  out  of  the  three  flight  conditions.  This 
resulted  in  a  controller  that  yielded  satisfactory  results  for 
the  0.6  mach,  sea  level,  and  0.9  mach,  30,000  ft.  flight  condi¬ 
tions.  This  universal  controller  consists  of  the  controller  used 
for  the  0.9  mach,  30,000  ft.  flight  condition.  Interestingly, 
this  corresponds  to  Bauschllcher' s  use  of  the  longitudinal 
controller  from  the  0.9  mach  flight  condition  as  his  universal 
longitudinal  controller. 

Not  having  a  universal  controller  for  all  three  flight 
conditions  may  not  be  a  detriment  in  this  case,  because  of  the 
unique  nature  of  the  0.15  mach,  sea  level  flight  condition.  This 
flight  condition  is  expected  to  be  different  from  the  other  two 
due  to  the  very  slow  airspeed  and  higher  angle  of  attack.  The 
controller  found  here  would  have  to  be  tested  with  other  more 
normal  flight  conditions  to  verify  its  use  as  a  universal 
controller . 

Testing .  Testing  the  universal  controller  consisted  of 
performing  the  longitudinal  maneuvers  for  the  0.6  mach,  sea  level 
flight  condition  using  the  controller  matrices  and  M  matrix  from 
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the  0.9  mach,  30,000  ft.  flight  condition.  The  same  maneuvers 
from  Chapter  IV  are  used  and  consist  of: 

Pitch  Pointing 

Direct  Lift 

Vertical  Translation 

Figures  44  through  52  show  the  results  of  these  maneuvers,  along 
with  the  commanded  inputs  and  control  surface  deflections.  Note 
that  the  title  for  each  plot  uses  the  abbreviation  "U"  to 
designate  universal  controller  testing.  Table  IX  summarizes  the 
results  of  these  plots. 

Results .  Generally,  the  results  of  the  responses  for  the 
universal  controller  are  fairly  good.  This  section  details  the 
results  for  each  maneuver. 

Figures  44  through  46  show  the  responses  for  the  1.75  degree 
pitch  pointing  command.  Comparing  these  plots  with  the  ones  for 
the  responses  of  the  aircraft  with  the  controller  designed  for 
this  flight  condition  (Figures  21  through  23),  the  worst  feature 
of  this  maneuver  is  the  decreased  tracking  of  the  angle  of  at¬ 
tack.  However,  there  is  less  overshoot  in  the  theta  response  but 
with  greater  error  after  about  five  seconds.  The  deviation  of 
the  flight  path  angle  is  greater  with  the  universal  controller 
than  without.  The  velocity  deviation  from  zero  is  greater  with 
the  universal  controller,  but  the  frequency  of  oscillation  is 
much  less,  and  the  damping  is  greater.  The  control  surface 
deflections  aren't  much  different  with  the  universal  controller, 
and  there  is  a  significant  decrease  in  both  the  required  magni¬ 
tude  of  thrust  and  its  rate  of  oscillation. 

Making  the  same  comparisons  of  the  universal  controller 
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1.00  2.00  4.00  0.00  0.00  10.00  12.00  14.00  10.00 

_ TIME.  SECONDS _ 

1.75  DEGREE  PITCH  POINTING  (0.6  MflCH.  0  FT) 

Figure  44.  1.75  Degree  Pitch  Pointing  With  Universal  Controll 

(0.6  Mach,  0  Ft.) 


COMMAND  INPUTS  ( THETA. U. ALPHA ) 


Figure  47 


1.8  G  Direct  Lift  With  Universal  Controller 
(0.6  Mach,  0  Ft.) 


U  (FT/SEC) 

r4.00  -3.00 


1.8  G  DIRECT  LIFT  (U.  0.6  MfiCH.  0  FT) 


1.8  G  DIRECT  LIFT  (U.  0.6  MflCH.  0  FT) 


Figure  48.  1.8  G  Direct  Lift  With  Universal  Controller 

(0.6  Mach ,  0  Ft  .  ) 
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1.8  0  DIRECT  LIFT  (U.  0.6  MflCH.  0  FT) 


Figure  49.  1.8  G  Direct  Lift  With  Universal  Controller 

(0.6  Mach,  0  Ft .  ) 


CONTROL  SURFACE  INPUTS 


TABLE  IX 

SUMMARY  OF  OUTPUT  RESPONSES  FOR 

UNIVERSAL  CONTROLLER 

Fit. 

Cond . 

Maneuver 

Command,  V 

Peak 

Time  To 

l 

Final  j 

Value 

Peak 

Value 

0.60 

Mach 

Pitch 

.5,  .0305,20,20 

.0354 

8.1  sec 

.0350 

0  ft 

• 

Pointing 

0,0, 0,0 

-.849 

2.1  sec 

.277 

.5, .0305,20,20 

.015 

10  sec 

.0149 

Direct 

10, .865,20,20 

.776 

10  sec 

.776 

Lift 

0,0, 0,0 

-3.93 

9.6  sec 

-3.83 

0,0, 0,0 

.0237  . 

1.9  sec 

.0168 

Vertical 

0,0, 0,0 

.00828 

2.4  sec 

.00597 

Transl . 

0,0, 0,0 

.403 

5  sec 

.203 

1, .0384,20,20 

.0189 

10  sec 

.0189 

The  peak  value  Is  the  greatest  deviation  from  the  Initial 
value  of  zero.  The  command  vector  Is  composed  of: 

~S  (rad)  “1 
V  ■  u  (ft/sec)i 
o<  (rad)  J 


TIME.  SECONDS 


0.8  G  VERTICAL  TRANSLATION  (U.  0.6  MACH.  0  FT) 

Figure  52.  0.8  G  Vertical  Translation  With  Universal  Controller 

(0.6  Mach ,  0  Ft . ) 
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response  (Figures  47  through  49)  with  the  responses  using  the 
controller  designed  for  this  flight  condition  (Figures  24  through 
26)  for  a  1.8  g  direct  lift  maneuver,  several  conclusions  can  be 
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drawn.  The  theta  response  with  the  universal  controller  is 
practically  Identical  to  the  dedicated  controller  response,  but 
once  again  the  velocity  response  has  greater  magnitude  of  devia¬ 
tion  but  with  less  frequency  of  oscillation.  The  angle  of  attack 
response  is  also  almost  identical  for  both  controllers.  The 
control  surface  deflections  along  with  the  thrust  required  are 
not  too  dissimilar.  Of  the  three  maneuvers,  this  is  the  one  with 
the  best  response  using  the  unviversal  controller. 

Comparing  the  responses  of  the  universal  controller  to  a  0.8 
g  vertical  translation  (Figures  50  through  52)  with  those  using 
the  dedicated  controller  (Figures  27  through  29)  yields  the  worst 
performance  for  the  universal  controller.  Even  though  the  theta 
response  has  improved,  the  angle  of  attack  response  has  signifi¬ 
cantly  degraded.  Once  again  the  forward  velocity  response  with 
the  universal  controller  exhibits  greater  deviation  from  the 
commanded  zero  level  but  has  greater  damping  and  lower  frequency. 
And  as  with  the  other  maneuvers  the  control  surface  deflections 
have  not  changed  much,  although  the  thrust,  in  addition  to  being 
smaller  in  magnitude,  also  has  a  marked  change  in  the  basic  shape 
of  the  time  response  plot. 

Conclusions .  Even  though  the  controller  from  the  0.9  mach 
flight  condition  proved  to  yield  the  best  results  of  any  of  the 
controllers  tested  in  the  search  for  a  universal  controller,  its 
performance  is  probably  not  adequate  enough  for  a  fighter 
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aircraft,  such  as  the  FPCC.  However,  there  is  some  degree  of 
robustness  to  the  universal  controller,  and  perhaps  this  type  of 
approach  would  yield  an  adequate  universal  controller  for  a 
transport  type  aircraft. 

Effects  of  Controller  Delay 

This  section  details  the  testing  of  the  effects  of 
introducing  a  one  sampling  period  time  delay  into  the  output  of 
the  PI  controller.  This  delay  represents  a  real  world  delay  that 
would  exist  in  the  output  of  a  digital  controller  using  a  sample 
and  hold  device  (such  as  an  analog  to  digital  converter)  to 
produce  the  input  error  signal.  This  is  implemented  in  the 
simulation  within  the  program  MULTI  by  simply  delaying  the  output 
of  the  controller  one  sampling  period,  i.e.  the  current  control 
input  to  the  actuators  (or  plant)  is  the  previously  calculated 
value,  while  the  most  recently  calculated  value  is  saved  for  the 
next  sampling  period  (which  is  the  "outer  loop"  within  the  MULTI 
simulation).  MULTI  also  has  the  capability  of  allowing  the 
analog  plant  (the  aircraft  in  this  thesis)  to  continue 
calculating  new  states  and  outputs  by  setting  the  calculation 
step  size  smaller  than  the  sampling  period.  Using  this  feature 
would  certainly  make  the  simulation  more  true  to  life  and  should 
be  considered  for  future  thesis  work. 


Results . 

Figures  53 

through  79  show 

the 

results  of 

performing  the 

longitudinal 

maneuvers  for  all 

of 

the  flight 

conditions  from  Chapter  IV  with  a  one  sampling  period  time  delay 
in  the  output  of  the  PI  controller.  Note  that  this  delay  is  only 
valid  for  the  case  of  a  digital  controller.  Table  X  summarize'. 


FORWARD  VELOC 
.60  -1.20  -0. 


THRUST  (LBS 


COMMAND  INPUTS  ( THETfi. U. ALPHA  ) 


FORHflRO  VELOCITY.  U  (FT/SEC)  _  I  CONTROL  SURFACE  INPUTS 


TABLE  X 


SUMMARY  OF  OUTPUT  RESPONSES  FOR  LONGITUDINAL  CONTROLLERS  WITH 

DELAY 


0.15  Mach 
0  ft. 


0.6  Mach 
0  ft. 


0.9  Mach 
30,000  ft. 


Maneuver 

C  ommand ,  V 

Peak 

Time  To 

Final 

Value 

Peak 

Value 

Pitch 

.5, .0305,20,20 

.0351 

1  sec 

.0310 

Pointing 

0,0, 0,0 

-.685 

0.7  sec 

-.0497 

.5, .0305,20,20 

.0336 

0.9  sec 

.0238 

Direct 

10, .976,20,20 

.950 

10  sec 

.950 

Lift 

0,0, 0,0 

-1.58 

0.9  sec 

-.651 

0,0, 0,0 

.0924 

3.4  sec 

.0141 

| 

Vertical 

0,0, 0,0 

.00244 

2.6  sec 

.000663  | 

Transl . 

0,0, 0,0 

.289 

1.2  sec 

-.000325i 

1  ,  .0977, 20, 20 

.0824 

10  sec 

.0824 

>wing  flight  conditions  reflect  delay  compensation)  ; 

Pitch 

.5, .0305,20,20 

.0571 

1.6  sec 

i 

.0416 

Pointing 

o 

o 

w 

o 

o 

-.370 

0.8  sec 

.0893 

.5, .0305,20,20 

.0248 

1  sec 

.0236  | 

Direct 

10, .865,20,20 

.774 

10  sec 

1 

.  774  i 

Lift 

0,0, 0,0 

-.583 

8.5  sec 

-.540  ! 

0,0, 0,0 

.0183 

1.5  sec 

.00723  j 

Vertical 

0 ,0,0,0 

.0329 

2  sec 

.0149 

Transl . 

0, 0, 0, 0 

.242 

4.1  sec 

.0362 

1, .0384,20,20 

.0299 

10  sec 

.0299 

Pi  tch 

.5, .0524,20,20 

.0635 

2.1  sec 

.0569  j 

Pointing 

0,0, 0,0 

-1.516 

1.9  sec 

-  110  j 

.5, .0524,20,20 

.0267 

1.2  sec 

.0153 

Direct 

10, .865,20,20 

.784 

10  sec 

.784 

Lift 

o,o,e,o 

-5.50 

9.9  sec 

-5.49 

0,0, 0,0 

.0459 

2.4  sec 

.0361 

Vertical 

0,0, 0,0 

-.005 

2.9  sec 

-.00407 

Transl . 

0,0, 0,0 

-.187 

4.9  sec 

-.0852 

1,-. 036, 20, 20 

-.0107 

10  sec 

-.0107 

The  delay  compensation  Is  described  in  text.  The  peak 
value  Is  the  greatest  deviation  from  the  initial  value  of 
zero.  The  command  vector  is  composed  of: 

Q  (rad) 

V  -  u  (ft/ sec ) 

<*  (rad) 
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TABLE  XI 

SUMMARY  OF  OUTPUT  RESPONSES  FOR  UNIVERSAL  CONTROLLER  WITH  DELAY 


Fit . 

Cond . 

Maneuver 

Command,  V 

Peak 

Time  To 

Final 

Value 

Peak 

Value 

0.60 
0  ft 

Mach 

• 

Pitch 

Pointing 

.5, .0305,20,20 
0, 0,0,0 
.5, .0305,20,20 

.0369 

-.938 

.0130 

2.1  sec 
2.3  sec 

1  sec 

.0358 

-.0398 

.010 

Direct 

Lift 

10, .865,20,20 
0,0, 0,0 

0,0, 0,0 

.781 

-5.60 

.0259 

10  sec 

10  sec 
2.1  sec 

.781 

-5.60 

.0212 

Vertical 
Transl • 

0,0, 0,0 

0,0, 0,0 

1, .0384,20,20 

.00906 

-.312 

.0126 

2.6  sec 

2  sec 

10  sec 

.00749 

.0586 

.0126 

Figures  shown  represent  compensated  results,  with 
compensation  used  explained  in  the  text .  The  peak  value 
is  defined  as  greatest  deviation  from  the  initial  zero 
value.  The  command  vector  is  composed  of: 

~Q  (rad) 

V  -  u  (ft/sec) 

<x  (rad) 


TIME.  SECONDS 


0.5  G  VERTICAL  TRANSLATION  (0.  0 . ! 5  MACH.  0  FT) 

Figure  61.  0.5  G  Vertical  Translation  With  Delay 

(0.15  Mach,  0  Ft.) 


RAO 


0.00 

2.00 

4.00 

0.00  0.00  10.00  12.00 

14.00 

10.00 

TIME.  SECONDS 

1.75  DEGREE  PITCH  POINTING  ID.  ALPHA=2.  0.6  MACH J 


1.75  OEGREE  PITCH  POINTING  (0.  ALPHA=2.  0.6  MACH) 

Figure  62.  1.75  Degree  Pitch  Pointing  With  Delay,  Control  Ratio 

Of  Two  (0.6  Mach,  0  Ft.) 
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THETfl  ( RflD 


THRUST  l LBS)  .  ..  ...  *10 


1.8  G  DIRECT  LIFT  (0.  RLPHR=2.  0.6  MRCH ) 


1.8  0  OIRECT  LIFT  (0.  RIPHR=2.  0.6  MRCH) 


Figure  66.  1.8  G  Direct  Lift  With  Delay,  Control  Ratio  Of 

(0.6  Mach,  0  Ft . ) 


1.8  DIRECT  LIFT  (D.  ALPHA=2.  0.6  MACH) 


vvs 


TOTAL  THRUST 


0.8  0  VERTICAL  TRANSLATION  (0.  ALPHA=2 .  0.6  MACH) 


Figure  71.  0.8  G  Vertical  Translation  With  Delay,  Control  Ratio 

Of  Two  (0.6  Mach,  0  Ft.) 
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THRUST  UBS)  «10'  I  CONTROL  SURFACE  INPUTS 


RAO 


LBS 


TIME.  SEC0NOS 


1  0  VERTICAL  TRANSLATION  (0.  ALPHA* 4.  0.9  HACH) 


Figure  78 


1  G  Vertical  Translation  With  Delay,  Control  Ratio 
Of  Four  (0.9  Mach,  30000  Ft.) 


these  results 


Starting  with  the  pitch  pointing  maneuver  at  the  0.15  mach 
flight  condition  (Figures  53  through  55),  the  effects  of  the 
delay  are  most  obvious  in  the  angle  of  attach  response.  Both 
theta  and  alpha  remain  essentially  unchanged  in  their  response 
during  the  first  1.5  seconds  when  compared  with  the  corresponding 
responses  for  no  delay  (Figures  12  through  14).  After  the 
initial  1.5  seconds  the  alpha  response  falls  off  considerably. 
The  forward  velocity,  control  surface  deflections,  and  total 
thrust  do  not  change  much  with  the  addition  of  the  delay. 

For  the  direct  lift  maneuver  at  0.15  mach  the  introduction 
of  the  delay  causes  almost  no  change  in  the  output  responses 
(Figures  56  through  58)  as  compared  to  the  same  responses  with  no 
delay  (Figures  15  through  17).  For  the  control  surfaces  the 
effect  of  the  delay  is  minimal,  as  is  the  effect  of  the  delay 
upon  the  total  thrust.  Interestingly,  this  maneuver  exhibits 
both  the  best  responses  with  and  without  the  delay. 

Incorporation  of  the  delay  into  the  vertical  translation 
maneuver  at  0.15  mach  (Figures  59  through  61)  does  degrade  the 
alpha  response  significantly,  although  the  theta  response  is 
essentially  unchanged  when  compared  to  the  responses  without  the 
delay  (Figures  18  through  20).  The  forward  velocity,  control 
surface  deflections,  and  total  thrust  responses  have  not  changed 
much  with  the  added  delay. 

Figures  62  through  63  show  the  results  of  adding  the  delay 
to  the  pitch  pointing  maneuver  for  the  0.6  mach,  sea  level  flight 
condition.  The  added  instability  from  the  delay  has  caused  the 
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theta  and  alpha  response  to  become  oscillatory,  and  the  forward 


velocity  response  has  become  unstable.  The  control  surface 


deflections  show  the  oscillatory  nature  and  the  total  thrust 


demonstrates  instability 


Obviously  these  results  are 


unacceptable  compared  to  those  without  any  delay  (Figures  21 


through  23). 


However,  if  the  ratio  of  proportional  to  integral  controller 


matrices  is  Increased  to  three  from  the  original  value  of  two. 


then  most  of  the  effect  of  the  delay  can  be  removed.  Figures  64 


through  65  show  the  responses  for  the  pitch  pointing  maneuver  at 


0.6  mach,  but  with  this  ratio  set  to  three.  Compared  with  the 


undelayed  responses  the  theta  and  alpha  plots  are  not  changed 


much,  nor  is  the  forward  velocity  response.  The  control  surface 


deflections  haven't  changed  much  either,  but  the  total  thrust 


response  has  changed.  The  combination  of  the  delay  and 


Increasing  the  proportional  to  Integral  control  ratio  has  yielded 


a  much  more  acceptable  thrust  response.  It's  important  to  point 


out  that  increasing  the  control  ratio  without  the  delay  does  not 


yield  satisfactory  results  for  all  of  the  maneuvers  at  this 


flight  condition. 


The  destabilizing  effect  of  the  delay  is  also  apparent  in 


the  responses  for  the  direct  lift  maneuver  at  0.6  mach  (Figures 


66  through  67).  The  theta  response  has  become  oscillatory,  with 


the  forward  velocity  once  again  showing  instability.  The  jet 


flaps  also  demonstrate  an  oscillatory  response,  with  the  total 


thrust  apparently  unstable.  These  results  would  have  to  be 


considered  unsatisfactory  when  compared  with  the  undelayed 


responses  (Figures  24  through  26). 
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Increasing  the  ratio  of  proportional  to  integral  control 
once  again  compensates  for  practically  all  of  the  deleterious 
effects  of  the  delay.  Figures  68  through  69  show  the  responses 
for  the  direct  lift  maneuver  at  0.6  mach  with  the  control  ratio 
increased  to  three  from  the  original  value  of  two.  The  theta, 
alpha,  and  forward  velocity  responses  are  nearly  identical  to  the 
undelayed  responses,  and  the  control  surface  deflections  are  very 
similar.  The  thrust  response  is  actually  Improved,  showing  less 
oscillatory  nature. 

The  vertical  translation  maneuver  at  0.6  mach  also  follows 
right  along  with  the  previous  results.  Figures  70  through  71 
show  the  responses  with  the  control  ratio  unchanged  but  with  the 
addition  of  the  delay.  The  theta,  alpha,  and  forward  velocity 
curves  show  unstable  oscillations  as  do  the  control  surface 
deflections  and  the  total  thrust  plots.  These  results  are  also 
unsatisfactory  when  compared  to  the  original  responses  (Figures 
27  through  29). 

Unfortunately,  increasing  the  control  ratio  to  three  from 
its  initial  value  of  two  does  not  have  quite  the  same  impact  as 
with  the  prior  maneuvers  (Figures  72  and  73).  The  theta  and 
alpha  responses  have  improved,  but  not  to  the  point  where  they 
coul.  be  called  similar  to  the  undelayed  responses.  They  are  not 
very  dissimilar  either,  but  the  effects  of  delay  are  not  as 
completely  compensated  for  as  before.  The  forward  velocity 
hasn't  changed  much,  nor  have  the  control  surface  deflections. 
But  once  again  the  thrust  response  has  actually  improved. 

Figures  74  through  79  show  the  responses  for  the 


longitudinal  maneuvers  at  the  0.9  teach ,  30,000  ft.  flight 
condition  with  attempted  compensation  thru  increasing  the  ratio 
of  proportional  to  integral  control.  These  results  can  be 
compared  to  the  undelayed  responses  for  this  flight  condition 
(Figures  30  through  39).  The  pitch  pointing  maneuver  is 
compensted  very  well,  with  an  actual  improvement  in  the  forward 
velocity  and  thrust  responses.  The  direct  lift  maneuver  is  also 
well  compensated,  but  without  any  obvious  improvement  in  the 
thrust  or  forward  velocity  responses.  The  vertical  translation 
also  is  well  compensated  without  any  marked  difference  in  any  of 
the  responses.  Unfortunately,  this  last  maneuver  is  one  of  the 
worst  responding  ones,  even  undelayed. 

Conclusions .  Due  to  time  constraints  in  preparing  the 
results  of  the  effects  of  and  compensating  for  time  delay,  all 
of  the  longitudinal  maneuvers  were  tested  instead  of  the  lateral 
maneuvers  at  the  0.15  mach  flight  condition.  This  was  done  so 
that  tentative  conclusions  could  be  drawn  regarding  the  effect  of 
delay  on  a  universal  controller. 

Obviously,  introducing  a  one  sampling  period  time  delay  into 
the  output  of  the  PI  controller  introduces  oscillations  and 
instabilities  into  the  time  responses.  The  effect  of  the  delay 
seems  to  increase  with  increasing  speed  and  altitude,  although 
this  is  a  tenuous  conclusion  at  best.  For  this  aircraft, 
compensating  for  this  time  delay  is  effectively  achieved  through 
Increasing  the  ratio  of  proportional  to  integral  control  in  the 
PI  controller  matrices.  Time  limitations  precluded  rigorous 
searching  for  other  methods  of  compensating  for  this  delay.  One 
of  the  most  interesting  results  of  compensating  for  the  delay  is 
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an  improvement  of  the  velocity  and/or  thrust  responses  for  some 
of  the  maneuvers.  Overall,  this  method  of  compensation  seems 
very  promising. 


Universal  Controller  Incorporating  Delay 

Extending  these  results  to  compensating  a  universal  type 
controller,  the  logical  conclusion  would  be  that  Increasing  the 
ratio  of  proportional  to  Integral  control  should  work,  with  the 
actual  Increase  probably  only  found  experimentally.  To  test  this 
hypothesis,  the  universal  controller  discussed  at  the  beginning 
of  this  chapter  is  tested  by  introducing  the  time  delay.  Figures 
80  through  85  show  the  time  responses  for  the  universal 
controller  using  the  longitudinal  maneuvers  at  the  0.6  mach 


flight  condition. 


Limited  experimentation  showed  that  a 


proportional  to  integral  control  ratio  of  four  was  close  to  the 
best  value,  with  Increasing  ratios  degrading  the  angle  of  attack 
response  and  decreasing  ratios  not  compensating  the  delay  enough. 
This  is  the  same  ratio  used  for  the  0.9  mach  flight  condition, 
further  reinforcing  the  choice  of  the  design  parameters  from  the 
0.9  mach  flight  condition  for  the  universal  controller.  Table  XI 
summarizes  the  results. 

Results .  Figures  80  through  81  show  the  time  responses  of 
the  universal  controller  to  the  1.75  degree  pitch  pointing  com¬ 
mand.  Comparing  these  to  the  results  for  the  undelayed  responses 
(Figures  44  through  46)  shows  no  major  changes  in  the  time 
responses.  The  theta  response  has  more  overshoot,  while  the 
angle  of  attack  shows  greater  error  at  the  end  of  ten  seconds. 
The  transient  response  of  alpha  is  almost  unchanged.  The  forward 
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Figure  81 


1.75  Degree  Pitch  Pointing  With  Delay,  Universal 
Controller  (0.6  Mach,  0  ft.) 


2.4  G  DIRECT  LIFT  (U.  0.  ALPHfl=4.  0.6  M ) 


2.4  G  OIRECT  LIFT  (U.  0.  ALPHfl=4.  0.6  M ) 

Figure  82.  1.8  G  Direct  Lift  With  Delay,  Universal  Controller 

(0.6  Mach,  0  ft.) 
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Figure  84 


0.8  G  Vertical  Translation  With  Delay,  Univers 
Controller  (0.6  Mach,  0  ft.) 
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velocity  shows  increased  damping  but  with  about  the  same  magni¬ 
tude  of  deviation  from  zero.  The  control  surface  deflection  for 
the  horizontal  canards  hasn't  changed  much,  while  the  jet  flaps 
are  not  deflected  as  much.  Total  thrust  requirements  are  reduced 
slightly . 

From  Figures  82  through  83  comparisons  can  be  made  with  the 
undelayed  responses  (Figures  47  through  49)  for  the  1.8  g  direct 
lift  command.  Once  again,  the  responses  have  changed  little  with 
the  introduction  of  the  delay  and  subsequent  compensation  for  the 
delay.  The  theta  and  alpha  responses  are  practically  identical, 
with  only  minor  changes  in  the  forward  velocity  response.  As 
with  the  pitch  pointing  command,  the  canard  deflection  has  not 
changed  much,  and  the  jet  flaps  are  deflected  less.  The  total 
thrust  responses  are  very  similar. 

The  final  maneuver,  the  0.8  g  vertical  translation  (Figures 
84  and  85),  yields  the  most  dissimilar  results.  The  alpha 
response  changes  quickly  after  the  first  second,  and  after  ten 
seconds  is  much  less  than  for  the  undelayed  response.  The  theta 
response  shows  more  overshoot,  but  otherwise  is  very  similar. 
Although  the  magnitude  of  the  forward  velocity  response  is  about 
the  same,  the  damping  is  much  greater  for  the  compensated  delayed 
plot.  Just  as  with  the  other  maneuvers  the  canard  deflection  is 
basically  unchanged  while  the  jet  flaps  are  deflected  less.  The 
thrust  requirements  are  only  slightly  reduced. 

Conclusions .  Once  again  the  overriding  conclusion  is  that 
if  the  maneuver  performs  well  without  the  delay,  then  the  delay 
can  be  effectively  compensated  for  by  increasing  the  ratio  of 
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proportional  to  Integral  control,  whereas  a  maneuver's 
performance  is  not  going  to  improve  with  the  inclusion  of  the 
delay.  As  with  other  results  from  this  thesis,  the  angle  of 
attack  response  seems  to  be  the  limiting  factor  in  trying  to 
improve  aircraft  response.  This  is  most  evident  in  the  pitch 
pointing  maneuver  where  Increasing  the  control  ratio  decreases 
the  alpha  response. 

Minimizing  Controller  Matrices  Elements 

This  section  describes  the  efforts  to  reduce  the  number  of 
elements  needed  in  the  controller  matrices  by  finding  common 
elements  to  be  set  equal  to  zero.  If  the  same  elements  in  either 
controller  matix  can  be  set  to  zero  for  all  maneuvers  at  a  given 
flight  condition,  then  for  that  flight  condition  fewer  memory 
locations  need  to  be  accessed,  increasing  computational  speed. 
This  is  especially  true  in  consideration  of  the  fact  that  a 
memory  access  is  one  of  the  slowest  microprocessor  instructions. 
The  greatest  benefit  would  occur  if  the  same  elements  could  be 
eliminated  for  all  flight  conditions  (or  perhaps  all  common 
flight  conditions).  Unfortunately,  the  later  case  was  not  true 
for  this  thesis  effort.  The  best  results  obtained  were  a 
reduction  of  two  elements  from  the  controller  matrices  for  the 
longitudinal  maneuvers  at  the  0.15  mach  flight  condition,  and  a 
reduction  of  one  element  for  the  lateral  maneuvers  at  the  same 
flight  condition.  This  means  that  a  combined  controller  for  the 
0.15  mach  flight  condition  would  need  three  fewer  elements.  For 
the  other  flight  conditions  no  elements  were  found  that  could  be 
set  to  zero  for  all  of  the  maneuvers.  In  fact,  for  the  0.9  mach 
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flight  condition,  there  were  no  elements  that  could  be  set  to 
zero  for  even  two  of  the  maneuvers  (which  was  possible  at  0.6 


mach)  .  Consequently,  the  following  discussion  covers  only  the 
controllers  for  the  0.15  mach  flight  condition. 

Lateral  Controller .  For  the  lateral  controller,  the  follo¬ 
wing  matrix  is  used  for  the  proportional  and  integral  controller 
matices  (which  are  the  same).  The  zero  element  represents  an 
element  that  could  be  removed  from  the  algorithm  used  within  the 


PI  processor. 
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Figures  86  through  87  show  the  responses  of  the  aircraft  to 
the  1.5  g  horizontal  translation  command  using  the  controller 
matrix  from  equation  (98).  Comparing  these  to  the  original 
responses  to  this  command  using  the  full  controller  matrices 
(Figures  5  through  6),  yields  some  interesting  observations. 
First,  the  side  velocity  response  is  almost  exactly  identical 
with  the  original  and  the  yaw  rate  response  has  changed  very 
little.  However,  the  roll  angle  response  has  improved 
considerably.  The  control  surface  deflections  show  very  similar 
characteristics  also.  This  means  that  the  maneuver  actually 
improved  by  eliminating  the  element  in  the  controller  matrices. 

From  Figures  88  through  89  similar  observations  are  made 
comparing  the  results  of  the  reduced  controller  matrices 
responses  to  the  1  g  flat  turn  to  the  original  unmodified 
responses  (Figures  7  through  8).  The  side  velocity  and  yaw 
rate  responses  have  changed  very  little  while  the  roll  angle 
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Figure  89.  1  G  Flat  Turn  With  Reduced  Controller  Matrices 

(0.15  Mach,  0  ft.) 
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3  DECREE  YflW  POINTING  (Z.  0.15  MACH) 


Figure  90.  3  Degree  Yaw  PoinCing  With  Reduced  Controller 

Matrices  (0.15  Mach,  0  ft.) 
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Figure  92 


1.75  Degree  Pitch  Pointing  With  Reduced  Controller 
Matrices  (0.15  Mach,  0  ft.) 
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Figure  94.  1.75  Degree  Pitch  Pointing  With  Reduced  Controller 

Matrices  (0.15  Mach,  0  ft.) 


Figure  95.  0.5  G  Direct  Lift  With  Reduced  Controller  Matrices 

(0.15  Mach,  0  ft.) 
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Figure  96.  0.5  G  Direct  Lift  With  Reduced  Controller  Matrices 

(0.15  Mach,  0  ft.) 
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Figure  97.  0.5  G  Direct  Lift  With  Reduced  Controller  Matrices 

(0.15  Mach,  0  ft.) 


0.5  G  VERTICAL  TRANSLATION  (Z.  0.15  MACH ) 


Figure  98.  0.5  G  Vertical  Translation  With  Reduced  Controller 

Matrices  (0.15  Mach,  0  ft.) 
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Figure  99 


0.5  G  Vertical  Translation  With  Reduced  Controller 
Matrices  (0.15  Mach,  0  ft.) 
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Figure  100.  0.5  G  Vertical  Translation  With  Reduced  Controller 

Matrices  (0.15  Mach,  0  ft.) 


response  has  been  improved  quite  dramatically.  Once  again  the 
control  surface  deflections  have  not  changed  except  for  their 
transient  response.  This  maneuver  has  also  been  improved  merely 
by  introducing  a  zero  element  in  the  controller  matrices. 

For  the  final  maneuver,  the  3  degree  yaw  pointing.  Figures 
90  through  91  show  that  very  little  changes  have  occured  in  any 
of  the  responses  when  compared  with  the  original  unmodified 
curves  (Figures  9  through  11).  In  fact,  the  roll  angle  response 
has  degraded  just  slightly. 

Since  the  element  set  to  zero  would  contribute  to  the 
aileron  deflection  depending  on  the  error  in  yaw  rate,  then  the 
aileron  deflection  would  be  based  more  upon  roll  angle  error. 
This  should  reduce  the  roll  angle  error  since  roll  angle  is 
mainly  governed  by  aileron  deflection.  The  above  results 
certainly  seem  a  logical  extension  of  this  analysis. 
Unfortunately,  this  author  found  no  aprlori  method  of  predicting 
which  elements  of  the  controller  matrices  could  be  set  to  zero 
based  upon  this  type  of  analysis. 

Longitudinal  Controller .  The  following  shows  which  elements 
of  the  controller  matrices  have  been  set  to  zero.  Zero  elements 
can  be  removed  from  the  algorithm  that  is  used  within  the  PI 
controller  to  calculate  the  control  inputs  to  the  actuators. 
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Figures  92  through  94  show  the  responses  of  the  reduced 


element  controller  to  the  1.75  degree  pitch  pointing  command. 
Comparing  these  results  with  the  original  responses  with  no 
modifications  (Figures  12  through  14),  a  couple  of  observations 
are  possible.  The  most  striking  is  that  although  the  responses 
of  theta,  alpha,  forward  velocity,  and  control  surface 
deflections  are  almost  identical  with  the  original  plots,  the 
total  thrust  is  not.  In  fact,  with  the  introduction  of  the  two 
zero  elements,  the  maximum  required  thrust  has  been  reduced  by 
about  a  factor  of  four.  Why  this  would  be  with  no  apparent 
changes  in  any  of  the  other  responses  seems  to  have  no  obvious 
logical  explanation. 

Figures  95  through  97  show  the  responses  for  the  reduced 
element  controller  to  the  0.5  g  direct  lift  command.  Comparing 
these  to  the  responses  for  the  original  unmodified  controller 
(Figures  15  through  17),  some  observations  can  be  made.  Once 
again  the  theta  and  alpha  response  are  essentially  unchanged,  but 
now  the  forward  velocity  shows  a  surprising  change.  Whereas  with 
the  original  responses  the  velocity  response  seemed  unstable 
after  about  six  seconds,  now  the  velocity  response  shows  no  such 
instability.  There  seems  to  be  no  obvious  logical  reason  for 
this.  The  control  surface  deflections  seem  to  have  changed  very 
little,  but  close  inspection  of  the  curves  from  Figure  96  shows 
that  they  also  differ  from  the  original  curves  after  about  six 
seconds.  The  most  obvious  difference  is  the  thrust  response, 
which  has  a  marked  lack  of  the  oscillatory  and  possible  unstable 
response  of  the  original  response.  This  maneuver  has  definitely 
been  improved  in  the  steady-sta;e  response  by  introducing  the  two 


zeros 


Figures  98  through  100  show  that  the  instabilities  present 
in  the  original  responses  to  the  0.5  g  vertical  translation 
command  (Figures  18  through  20)  have  not  disappeared  with  the  two 
zero  elements  in  the  controller  matrices.  The  theta,  alpha,  and 
forward  velocity  responses  haven't  changed  very  much,  although 
the  velocity  response  shows  slightly  smaller  peak  values.  The 
control  surface  deflections  show  little  change,  but  the  thrust 
response  has  changed  significantly.  Unfortunately,  this  maneuver 
shows  no  signs  of  improved  response  with  the  two  zero  elements, 
although  there  appears  to  be  no  degradation  of  the  responses. 

The  two  zero  elements  reduce  the  contributions  to  the 
horizontal  canard  deflection  to  just  the  pitch  angle  error. 
Following  the  train  of  thought  from  the  lateral  case,  it  would  be 
expected  that  if  anything  would  be  improved  it  would  be  the  theta 
response,  and  yet  this  did  not  happen.  However,  the  one  zero 
that  corresponds  to  velocity  error  contributions  to  canard 
deflections  is  an  obvious  candidate  for  elimination.  But  the 
other  zero,  the  angle  of  attack  error  effect  on  canard  deflection 
is  not  such  an  obvious  choice.  Also,  other  seemingly  obvious 
choices  for  zero  elements  did  not  work  out. 

Conclusions .  Even  though  there  seems  to  be  no  apriori  way 
to  pick  the  zero  elements  in  the  controller  matrices,  the 
possible  choices  will  always  be  finite  and  can  be  tested  using 
intelligent  software.  The  expected  benefits  from  doing  this 
include  possibly  improved  responses  and  the  reductions  of  memory 
accesses  needed  in  the  controller  algorithm.  These  factors  would 
definitely  recommend  this  technique  for  "polishing"  the  design. 


CHAPTER  VI 


RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Introduction 

This  thesis  incorporates  Professor  Porter's  multivariable 
design  method  into  the  design  of  flight  controllers  for  three 
flight  conditions  for  the  FPCC  aircraft.  The  first  flight 
condition  represents  a  possible  STOL-like  operation,  while  the 
other  two  are  more  representative  of  normal  operation  and  were 
also  studied  by  Jon  Bauschlicher  in  a  previous  thesis  (Ref  2). 
One  of  the  purposes  of  this  thesis  was  to  design  longitudinal 
controllers  for  these  two  flight  conditions  using  a  different 
longitudinal  model  than  Bauschlicher  so  as  to  make  comparisons. 

For  the  first  flight  condition,  0.15  mach  at  sea  level,  a 
complete  controller  is  designed  by  decoupling  the  aircraft  model 
into  a  longitudinal  and  a  lateral  model.  The  other  two  flight 
conditions,  0.6  mach  at  sea  level,  and  0.9  mach  at  30,000  ft., 
are  decoupled  into  a  lateral  and  longitudinal  model,  but  only  the 
longitudinal  design  is  accomplished  because  Bauschlicher  had 
already  designed  controllers  for  the  lateral  model.  An  attempt 
to  find  a  "universal"  controller  capable  of  performing  all  the 
maneuvers  for  all  of  the  flight  conditions  failed,  but  the 
controller  for  the  0.9  mach  flight  condition  is  found  to  yield 
marginal  results  at  the  0.6  mach  flight  condition. 

A  digital  implementation  of  the  PI  controller  would 
introduce  a  delay  in  the  output  of  the  controller,  and  the 
effects  of  this  delay  are  studied  for  all  of  the  maneuvers  at  all 
of  the  flight  conditions.  Finally,  an  attempt  to  find  common 


elements  of  the  controller  matrices  that  can  be  eliminated  Is 


partially  successful. 

Additionally,  the  design  tool  MULTI,  a  program  based  upon 
Professor  Porter's  design  ‘method  and  capable  of  simulating 
aircraft  response  to  commanded  maneuvers,  is  upgraded  with  two 
more  options  added. 

Results 

Figures  5  through  11  show  the  time  responses  for  the  outputs 
and  the  flight  control  surface  deflections  for  the  0.15  mach 
flight  condition  lateral  controller.  The  figures  of  merit  and 
the  commanded  inputs  for  each  maneuver  are  listed  in  Table  VII. 
One  controller  was  used  for  all  of  the  maneuvers  yielding  results 
that  were  not  much  different  than  if  an  optimum  controller  had 
been  used  for  each  maneuver.  Of  the  three  maneuvers,  the  yaw 
pointing  showed  the  best  results,  with  the  horizontal  translation 
fairly  acceptable.  ,  The  flat  turn  responses  would  have  to  be 
considered  unacceptable. 

Figures  12  through  20  show  the  time  responses  for  the 
outputs  and  surface  deflections  for  the  three  longitudinal 
maneuvers  at  0.15  mach,  sea  level.  Table  VIII  shows  the  figures 
of  merit  and  commanded  inputs  for  the  maneuvers.  Once  again  one 
controller  was  found  for  all  three  maneuvers,  and  once  again  the 
results  are  mixed.  The  pitch  pointing  maneuver  is  classified 
fair,  with  the  direct  lift  considered  very  good.  But  the 
vertical  translation  would  have  to  be  termed  poor. 


Figures  21  through  39  show  the  same  curves  for  the  other  two 
flight  conditions,  0.6  mach  at  sea  level,  and  0.9  mach  at  30,000 
feet.  Table  VIII  also  lists  the  figures  of  merit  and  commanded 


inputs  for  these  longitudinal  maneuvers.  At  0.6  mach,  the  direct 


lift  maneuver  shows  the  best  results,  which  would  have  to 
considered  very  good.  But  the  pitch  pointing  maneuver  would  be 
called  poor,  at  best,  and  the  vertical  translation  is  certainly 
unacceptable.  At  0.9  mach  the  results  are  about  the  same,  with 
the  direct  lift  showing  very  good  results  but  the  vertical 
translation  and  pitch  pointing  showing  unacceptable  results. 

The  general  conclusion  regarding  the  comparison  with 
Bauschlicher's  results  is  that  the  longitudinal  model  used  in 
this  thesis  did  not  perform  nearly  as  well  as  the  one 
Bauschlicher  used. 

An  attempt  to  find  a  universal  controller  failed,  and  the 
attempt  to  find  a  longitudinal  controller  that  would  work  for  two 
of  the  three  flight  conditions  yielded  marginal  results.  Figures 
44  through  52  show  the  responses  for  the  0.6  mach  flight 
condition  using  the  0.9  mach  controller.  The  direct  lift  showed 
very  good  results,  but  the  other  two  would  have  to  be  considered 
very  poor.  Table  IX  summarizes  the  results. 

All  of  the  longitudinal  maneuvers  for  all  three  flight 
conditions  were  repeated  with  a  one  sampling  period  time  delay 
added  to  the  output  of  the  PI  controller.  Figures  53  through  61 
show  these  results  for  the  0.15  mach  flight  condition.  As  can  be 
seen  from  all  of  the  maneuvers,  the  addition  of  the  delay 
introduces  instability  into  the  responses. 

Figures  62  through  73  show  the  effects  of  delay  on  the 
longitudinal  maneuvers  for  the  0.6  mach  flight  condition  both 
with  and  without  compensation.  The  compensation  is  an  Increase 
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in  the  ratio  of  proportional  to  integral  control  within  the  PI 
controller.  The  uncompensated  curves  show  quite  dramatically  the 
unstabilizing  effect  of  the  delay.  From  the  compensated  plots 
the  conclusion  is  that  the  delay  can  be  effectively  cancelled 
out . 

From  Figures  74  through  79  the  conclusion  is  that  the  delay 
can  also  be  compensated  for  in  the  longitudinal  maneuvers  at  the 
0.9  mach  flight  condition.  Table  X  summarizes  these  results  for 
all  of  the  longitudinal  maneuvers  at  the  three  flight  conditions. 

Extending  this  testing  to  the  semi-universal  controller 
yields  the  results  shown  in  Figures  80  through  85.  Once  again 
the  conclusion  is  that  the  increase  of  the  proportional  to 
integral  control  ratio  effectively  compensates  for  the  addition 
of  the  delay.  Table  XI  summarizes  these  results. 

The  final  testing  performed  is  an  attempt  to  reduce  the 
number  of  gains  needed  in  the  controller  matrices. 
Unfortunately,  there  were  no  common  gains  found  that  could  be 
eliminated  from  all  of  the  controller  matrices  for  all  three 
flight  conditions.  In  fact,  there  were  no  gains  that  could  be 
eliminated  from  the  matrices  for  all  three  maneuvers  at  either 
the  0.6  or  0.9  mach  flight  condition.  There  were  gains  that 
could  be  removed  from  both  the  longitudinal  and  lateral 
controllers  for  the  0.15  mach  flight  condition.  Figures  86 
through  100  show  that  these  reduced  gain  controllers  show  no 
degradation  in  response  and  even  sometimes  improve  the  response. 

Conclusions 


Summarizing  the  above  results,  the  conclusions  drawn  are 


these : 


1.  Removing  the  maneuver  flaps  from  the  longitudinal  model 
definitely  degrades  the  responses  to  the  commanded  maneuvers. 

2.  Adding  a  delay  to  the  output  of  the  PI  controller  adds 
instability  to  the  responses. 

3.  The  effects  of  the  delay  can  be  effectively  compensated 
for  by  increasing  the  ratio  of  proportional  to  integral  control. 

4.  Setting  certain  elements  of  the  controller  matrices  to 
zero  can  yield  improved  responses. 

In  addition  to  these  conclusions,  there  are  two  important 
results  worth  summarizing: 

1.  A  universal  controller  that  yielded  acceptable  results 
compared  to  the  Individual  controller  responses  was  not  found. 

2.  The  aircraft  responses  were  typically  limited  by  the 
angle  of  attack  response. 

Recommendations 

This  section  describes  the  recommendations  this  author  has 
for  future  efforts  involving  both  the  design  method  and  the 
design  tool  MULTI. 

Design  Method.  The  following  is  a  list  of  the  suggested 
areas  for  future  work  on  the  Porter  method  of  controller  design: 

1.  Investigate  the  transformation  from  the  analog  to  the 
discete  PI  controller  more  closely.  The  current  equations  for 
the  discete  controller  implement  the  integration  via  a  first 
order  rectangular  integration  approximation.  Perhaps  reworking 
these  equations  to  yield  more  accurate  integration  approximations 
would  be  valuable. 
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2.  More  information  about  the  relationship  between  the 
transmission  zeros  and  the  types  of  expected  results  is  needed. 
Although  Professor  Porter's  theory  implies  that  the  regular 
design  should  not  have  any  limitations  in  the  responses  due  to 
the  transmission  zeros,  this  author  found  that  in  general  an 
irregular  design  yielded  better  results.  If  irregular  designs 
are  (and  so  far  they  do  seem  to  be)  going  to  be  more  prevelant  in 
flight  control  design,  then  insight  into  how  transmission  zeros 
affect  responses  would  certainly  be  helpful.  An  example  of  the 
type  of  questions  that  could  be  answered  is:  How  do  transmission 
zeros  in  the  right  half  plane  or  at  the  origin  affect  controller 
design  and  expected  responses? 

3.  Since  irregular  designs  have  been  the  most  prevelant  to 
date,  most  of  the  design  work  has  centered  around  the  use  of 
angles  instead  of  rates  as  outputs.  This  author  feels  more 
designs  incorporating  rates  (and  therefore  typically  regular 
designs)  should  be  attempted.  Reasons  for  this  suggestion 
include  the  two  considerations  that  current  pilot  controls 
command  rates,  and  that  regular  designs  have  no  inherent 
limitation  upon  the  output  responses  due  to  the  slow  modes  (as  do 
irregular  designs).  Another  factor  is  that  for  an  initial 
aircraft  attitude  that  is  not  straight  and  level,  the  use  of 
angles  is  unappropriate. 

Design  Tool ,  MULTI .  The  following  is  a  list  of  suggestions 
for  improvements  to  the  design  tool  MULTI: 

1.  In  conjuction  with  the  suggestion  for  improving  the 
integration  approximation  in  the  design  method,  implement  an 
improved  integration  routine  within  the  simulation  in  MULTI. 
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Appendix  B  lists  the  program  code  for  one  such  improvement. 

2.  Consideration  should  be  given  to  implementing  the  exact 
discrete  differential  equations  instead  of  using  the  current 
analog  integration  approximations,  accomplished  with  the  system 
library  ODE  (ordinary  differential  equation  solver).  This  would 
not  be  a  formidable  task  because  the  system  marices  are  all  time 
invariant  (constant  coefficients).  Students  receive  the 
background  necessary  for  this  type  of  implementation  in  courses 
such  as  EE  7.12  and  EE  6.44.  If  the  simulation  were  written  in 
this  fashion  and  then  segmented  as  a  separate  overlay  it  would  be 
useful  as  a  general  purpose  discrete  simulation  of  analog 
systems . 

3.  The  calculation  of  actuator  (or  control  input)  rates 
should  be  added  to  Option  #28  in  MULTI,  along  with  a  toggle 
switch  that  the  user  can  use  to  specify  either  modified  spline  or 
linear  curve  fitting  on  the  CALCOMP  plots.  These  two  features 
used  together  would  allow  the  user  to  decide  when  to  use  linear 
curve  fitti-ig  based  on  curve  slopes.  The  program  code  to 
implement  these  two  features  is  included  in  Appendix  C. 

4.  MULTI  needs  to  be  completely  rewritten  in  a  language 
that  has  both  dynamic  memory  allocation  and  more  structure  than 
FORTRAN.  The  dynamic  memory  allocation  should  allow  for  larger 
numbers  of  states,  in  addition  to  being  more  efficient  in 
general.  Using  a  language  that  causes  program  code  to  have  more 
structure  should  help  prevent  what  has  happened  to  the  current 
version  of  MULTI,  which  is  a  proliferation  of  GO  TO  statements. 
Because  of  this  unstructured  nature  of  FORTRAN,  MULTI  has  evolved 


into  a  program  that  is  very  hard  to  understand  and  maintain. 


5.  The  addition  of  sensors  and  actuators  to  MULTI  (Options 
#4  and  #5)  has  caused  severe  problems.  Some  of  these  problems 
were  corrected  by  this  author,  but  a  couple  still  remain.  A 
quick  summary  of  the  known  remaining  problems  are  (all  are  within 
the  simulation): 

a.  Formation  of  the  y  measurement  vector  YM  is  no 
longer  correct  and  should  be  eliminated. 

b.  The  use  of  YM  to  form  the  error  vector  E  is  also  no 
longer  correct  and  the  YM  should  be  replaced  with  Q,  which  is  the 
output  vector  incorporating  sensors  (if  present).  This  is 
currently  causing  erroneous  data  for  a  regular  design  that 
Includes  sensors. 

c.  The  call  to  subroutine  YOUT  results  in  the  actual 
values  of  the  outputs  to  be  plotted  instead  of  the  sensed  values 
if  there  are  sensor?.  However,  this  is  not  clear  to  either  the 
programmer  or  the  user.  Addtionally,  this  means  that  there  is  no 
current  way  to  plot  the  sensed  values  of  the  outputs.  Something 
should  be  done,  even  if  it  is  just  improved  documentation  of  the 
program  code. 

6.  MULTI  should  be  changed  to  allow  the  user  to  specify 
both  the  number  and  placement  of  the  columns  in  the  M  matrix. 
This  would  permit  the  designer  to  chose  which  state  derivatives 
should  be  measured  (or  estimated).  Additionally,  the  provision 
for  an  M  matrix  should  be  added  to  the  options  for  a  regular 


design. 
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JUSTIFICATION  FOR  ELIMINATION  OF  JET  FLAPS  FROM  LATERAL  MODEL 


Introduction 

This  appendix  derives  the  relationship  between  the  jet  flaps 
and  the  derivatives  of  the  lateral  states  p  and  r,  showing  that 
the  jet  flaps  can  be  removed  from  the  lateral  model  used  In  this 
thesis.  This  relationship  Is  expressed  by  two  terms  In  the  B 
matrix  from  the  system  equation 
x  ■  Ax  +  Bu 

with  x  being  the  state  vector.  The  two  terms  are  In  the  first 
and  second  row,  sixth  column  of  the  original  12  state  system  B 
matrix  given  In  equation  (48)  from  Chapter  III.  From  equation 
(48)  these  two  numbers  are  given  as: 

.00939 

-.000238 

This  appendix  demonstrates  that  these  numbers  should  both  be 
zero,  removing  the  jet  flaps  from  the  lateral  partition  of  the  B 
matrix . 


,v.' 


Derivation  (Ref  4) 

Starting  with  the  rotation  versus  moment  equations 
0 

p  ■  Lg/Ixx  +  (non  linear  terms) 


r  -  Ng/Izz  +  (non  linear  terms) 


(A-l) 

(A-2) 


and  substituting 

L  ■  L  cos<*  -  N-slno^  +  Fy  A  z  +  FzgAy  -  Y 


^  ‘ 

Nfi  -  Lssin«  +  Nscoso<  +  Fxflay  +  FygAx  -  Xey  ^ 

^ 


(A-3) 


Klj  I 


(  A-4  ) 
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where  and  represent  the  aerodynamic  moments  in  stability 

axes.  Pulling  out  the  only  terms  that  are  dependent  upon  jet 
flap  setting: 

Fy8&z  -  0 
FzgA y  -  0 
Fxg/iy  -  0 
Fy0Ax  -  0 

since  &x.  Ay,  and  Az ,  the  components  of  the  center  of  gravity 
offset  from  the  moment  reference  point  are  all  set  to  zero  in  the 
FPCC  simulation  program.  Continuing  with  the  terms  that  are 
dependent  upon  jet  flap  setting,  the  direct  propulsion  force 

.  r,  CM  Sj  CosXAy  -  CDJ,  *4  f  VV*  5 

^  *  *lr  6  tojSj  cos  -  CDlj  >/3 S 

Yey  *  f,  CM^j  i 

Yeyi  —  ~V  fA  C’J  fj 

2e«y  *  ’*lf  fx 


Assuming  the  following  equalities  (which  are  true  in  the  FPCC 
simulation  program) 


with 


'&(yq  z 


I 


'  *  x  distance  from  center  of  gravity  to  engine  one 

XenyA  -  z  distance  from  center  of  gravity  to  engine  two 

ye^i  -  y  distance  from  center  of  gravity  to  engine  one 

-  y  distance  from  center  of  gravity  to  engine  two 


y  distance  from  center  of  gravity  to  engine  two 


•  z  distance  from  center  of  gravity  to  engine  one 
“  z  distance  from  center  of  gravity  to  engine  two 

*  yaw  offset  of  engine  one  with  respect  to  aircraft 

'Vci/ujA  "  yaw  offset  of  engine  two  with  respect  to  aircraft 

■  thrust  adjustment  factor  (0.98) 

■  scaler  total  velocity  (square  root  of  the  sum  of  the 
squares  of  the  body  axis  velocities) 

f*  ■  air  density 

S  -  planform  area 

Substituting  these  terms  back  into  equations  (A-3)  and  (A-4) 
results  in  neither  Ls  nor  Ng  having  any  dependency  on  the  jet 
flaps.  This  means  that  the  jet  flap  deflection  (considering  the 
above  stated  symmetry  and  engines  of  equal  thrust)  has  no  effect 
on  yaw  or  roll  rates,  a  logical  conclusion. 

Since  equations  (A-3)  and  (A-4)  substituted  into  (A-l)  and 
(A-2)  show  that  the  derivatives  of  p  and  r  have  no  dependency  on 
the  jet  flap  deflection,  then  these  corresponding  entries  in  the 
system  B  matrix  are  set  to  zero  for  the  purposes  of  this  thesis. 
Since  the  B  matrix  from  equation  (48)  shows  these  numbers 
Incorrectly  as  being  non  zero,  then  there  is  an  error  in  the  FPCC 
simulation  program.  Unfortunately,  this  means  that  any  other 
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APPENDIX  B 


IMPLEMENTING  A  FIRST  FORWARD  DIFFERENCE  APPROXIMATION  TO  THE 
I/S  ROUTINE  WITHIN  THE  PI  CONTROLLER  IN  MULTI 

Introduction 

The  simulation  program  MULTI  currently  uses  the  first 
backward  difference  approximation  for  the  Integration  routine 
within  the  proportional  plus  integral  (PI)  controller.  However, 
due  to  certain  results  In  this  thesis  (see  Figures  16,  19,  and 

20)  there  seemed  to  be  a  problem  with  the  Integration  within  the 
PI  controller.  The  referenced  plots  demonstrate  a  feature  that 
appears  to  be  the  result  of  windup.  Implying  that  the  Integration 
needed  improvement.  One  simple  way  of  improving  a  digital 
Integration  is  to  use  the  first  forward  difference  approxima¬ 
tion,  which  was  tried  with  the  routine  within  the  PI  controller 
with  good  results. 

Implementation 

The  current  integration  within  the  PI  controller  in  MULTI  is 
done  with  the  first  backward  difference  approximation  as 

expressed  in  the  line  of  FORTRAN  code  as: 

Z(I)-Z(I)+SAMPT*E(I)  (B-l) 

which  is  the  representation  of: 

Z[(k+l)T]-Z(kT)+TeI(k-l)T]  (B-2) 

Since  the  line  of  code  (B-l)  appears  after  the  calculation  of  the 
control  input  which  is  based  upon  Z(kT)  in 

u(kT)-{Koe(kT)+K,  Z(kT)}/T  (B-3) 

the  error  used  in  equation  (B-2)  is  actually  the  previous  error 


and  thus  Che  e[(k-l)T]  shown  in  equation  (B-2).  There  is  some 
question  as  to  whether  equation  (B-2)  is  actually  equivalent  to 
Professor  Porter's  expression  (Ref  xx) : 

Z  f (k+l)T]-Z(kT)+T[v(kT)-y(kT)]  (B-4) 

Close  examination  of  expression  (B-4)  shows  that  it  is  equivalent 
to  expression  (B-2)  if  the  next  iteration  of  kT  starts  at  the 
sample  and  hold  device  of  Figure  1.  However,  this  would  seem  to 
conflict  with  expression  (B-3)  where  the  error  used  to  generate 
the  control  input  is  in  the  same  iteration  as  the  Z.  One  way  to 
resolve  this  problem  would  be  to  say  that  the  kT  iterations  start 
with  the  sample  and  hold  device  and  then  to  define: 

Z[(k+l)T]-Z(kT)+Te(kT) 

which  1 8  a  first  forward  difference  approximation. 

What  was  tested  first  though,  was  an  average  of  the  first 
forward  and  first  backward  differences,  which  is  still  a  rectangular 
approximation,  thus  keeping  the  "spirit’’  of  Professor  Porter's 
work.  This  averaged  first  difference  can  be  expressed  as: 

Z[(k+l)TJ-Z(kT)+T{e(kT)+eI(k-l)TJ}/2  (B-5) 

The  FORTRAN  code  Implementation  of  expression  (B-5)  appeared  in 
MULTI  as: 

Z ( I )«Z ( I )  +SAMPT* ( E ( I ) +PR EVE (I))/2.0  (B-6) 

Note  that  in  expressions  (B-l)  and  (B-6)  the  running  variable  I 
is  used  to  position  within  the  vectors  involved,  not  as  a  kT 
iteration  counter,  and  that  PREVE  is  the  previous  error. 

With  the  code  (B-6)  substituted  into  the  PI  controller 
routine  the  plots  for  the  longitudinal  maneuvers  at  the  0.15  mach 
flight  condition  were  regenerated.  The  only  noticable  changes 
were  in  the  plots  for  the  forward  velocity  with  the  direct  lift 


command,  and  in  the  forward  velocity  and  total  thrust  curves  with 
the  vertical  translation  command.  Comparing  the  new  curves 
(Figures  B-l,  B-2 ,  and  B-3)  with  the  originals  (Figures  16,  19, 
and  20),  the  effect  of  using  the  new  integration  approximation  is 
to  eliminate  the  instabilities  near  the  end  of  the  simulation. 


After  testing  the  averaged  first  difference  approximation,  a 
first  forward  difference  approximation  was  implemented  by 


changing  (B-6)  to 

Z ( I ) -Z ( I ) +SAMPT*E ( I ) 


(  B-7  ) 


This  line  of  code  was  also  moved  before  the  calculation  of  the 
control  input  so  that  it  does  implement  the  first  forward 
difference  approximation.  Limited  testing  indicates  that  this 
change  also  removes  the  instabilities  due  to  using  the  first 
backward  difference  approximation. 


0.00  t.OO  4.00  S.00  0.00  10.00  12.00  14.00  10.00 

TIME.  SECONDS 


FIGURE  B-l 

Figure  B-l.  Forward  Velocity  Response  to  0.5  G  Direct  Lift 
Command  (0.15  Mach,  0  ft.) 
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APPENDIX  C 

ADDITIONS  TO  DESIGN  PROGRAM  MULTI 


Introduction 

This  appendix  briefly  describes  the  two  additions  made  to 
the  design  program  MULTI  by  this  author.  The  first  addition 
allows  the  use  of  the  local  file  names  for  the  system,  design, 
and  simulation  data  files  as  command  line  arguments.  Within 
MULTI,  these  command  line  arguments  are  then  used  in  an  Option 
#199  which  performs  the  same  functions  as  Options  #9,  #19,  and 
#29.  The  second  addition  is  Option  #28,  which  calculates  and 
displays  several  figures  of  merit  from  the  latest  simulation.  At 
the  end  of  this  appendix  is  a  short  section  on  the  addition  of 
program  code  that  would  calculate  and  display  maximum  actuator 
rates.  The  final  section  is  a  description  of  another  option  that 
would  allow  the  user  to  toggle  the  curve  fitting  for  the  CALCOMP 
plots  from  modified  spline  to  linear  and  back. 

Command  Line  Arguments  and  Option  #199 

The  concept  behind  this  addition  is  to  allow  the  user  to 
specify  the  local  file  names  of  the  three  files  containing  the 
system,  design,  and  simulation  parameters  as  command  line 
arguments  to  be  used  with  Option  #199.  For  example,  if  the  user 
had  these  three  files  named  MEMO,  MEM10,  and  MEM20,  then  the 
command  (in  response  to  Cyber's  Command  query) 

COMMAND>  MULTI ( , , MEMO , MEMl 0 , MEM20) 


will  start  up  MULTI  execution  in  the  normal  manner,  but  will 
allow  the  user  to  type  (in  response  to  the  Option  query) 


OPTION,  PLEASE  >  #199 

invoking  Option  #199.  The  user  must  type  in  all  three  local  file 
names  to  use  Option  #199.  If  not,  then  a  FORTRAN  file  read  error 
upon  attempting  Option  #199  will  cause  the  user  to  be  "bombed 
out"  of  MULTI  and  back  into  the  Cyber  command  mode.  Of  course, 
the  three  local  files  must  contain  the  data  in  the  format 
necessary  to  satisfy  MULTI's  requirements.  This  currently  means 
that  the  files  must  be  in  the  format  specified  in  the  MULTI  Users 
Manual  for  Options  #9,  #19,  and  #29  ’"his  manual  can  be  obtained 
from  the  Electrical  Engineering  De;  .tment,  Air  Force  Institute 
of  Technology,  Wright  Patterson  .  .  ’’orce  Base,  Ohio,  45433. 

Invoking  Option  #199  actually  causes  modified  versions  of 
Options  #9,  #19,  and  #29  to  be  sequentially  executed.  This  means 
that  the  user  can  use  Option  #199  whenever  the  other  three  would 
be  used.  The  following  shows  what  the  user  will  see  after 
specifying  Option  #199: 

OPTION,  PLEASE  >  #199 

DATA  COPY  COMPLETE  FOR  OPTIONS 
#2,  #3,  #4 

DATA  COPY  COMPLETE  FOR  OPTIONS 
#11,  #12,  #13,  #14,  #16,  #18 

DATA  COPY  COMPLETE  FOR  OPTIONS 
#21,  #22,  #23,  #24,  #25,  #27 

OPTION,  PLEASE  >  # 

At  this  point  the  user  may  continue  as  if  Options  #9,  #19,  and 

#29  had  just  been  completed. 

Local  file  names  as  command  line  arguments  is  not  currently 
supported  in  CDC  versions  of  FORTRAN  4.  If  MULTI  were  recompiled 
in  FORTRAN  4,  changing  the  program  declaration  line  in  the  main 


PROGRAM  EXEC  ( PLOT , TAPE  9-PLOT ) 

would  be  sufficient  (not  considering  other  syntax  problems)  to 
disable  the  command  line  argument  and  Option  #199  feature.  The 
code  for  Implementing  Option  #199  would  remain,  so  that  changing 
the  program  declaration  back  to: 

PROGRAM  EXEC  ( PLOT , TAPE  9-PLOT , TAPE 1 0 , TAPE 35 , TAPE 40 ) 
would  then  reactivate  command  line  arguments  and  Option  #199 
capability . 

Option  #28 ,  Figures  of  Merit 

The  second  addition  made  to  MULTI  by  this  author  is  the 
calculation  and  display  of  several  figures  of  merit  using  Option 
#28.  Specifically,  this  option  calculates  the  peak  value,  time 
to  peak  value,  minimum  value,  time  to  minimum  value,  final  value, 
and  settling  time.  All  of  the  calculations  are  performed  upon 
the  packed  data  used  for  the  plots.  Additionally,  Option  #28 
checks  to  see  if  the  simulation  has  been  run  prior  to  its 
execution,  and  if  not,  then  informs  the  user  and  returns  to  the 
option  query.  Assuming  the  simulation  has  been  performed,  then 
using  Option  #28  would  produce: 

OPTION,  PLEASE  >  #28 

THIS  OPTION  CALCULATES  THE  FIGURES  OF  MERIT 
CONTINGENT  UPON  COMPLETION  OF  SIMULATION 

HOW  MANY  SEQUENTIAL  OUTPUTS  (STARTING  WITH  Yl) 

DO  YOU  WISH  FIGURES  OF  MERIT  FOR?  >1 
Yl  PEAK-42.98407352055 

Yl  TIME  TO  PEAK-9.999999999995 
Yl  MINIMUM-0. 

Yl  TIME  TO  MINIMUM-0. 

DO  YOU  WISH  TO  CONTINUE  WITH  SETTLING  TIME 
CALCULATIONS?  (1  FOR  YES/0  FOR  N0)>1 


FINAL  VALUE  OF  Yl-42 . 98407352055 


DO  YOU  WISH  TO  USE  DEFAULT  VALUE  OF  WITHIN  2Z 

OF  THIS  VALUE  FOR  CALCULATION  OF  SETTLING  TIME? 

(1  FOR  YES/0  FOR  N0)>1 

SETTLING  TIME  FOR  YI-6 . 499999999999 
The  user  in  this  example  only  specified  one  output,  but  up  to  as 
■any  as  there  actually  are  can  be  specified,  with  the  values 
given  above  being  repeated  for  each  output.  Additionally,  the 
user  does  not  have  to  continue  with  settling  time  calculations, 
in  which  case  the  user  is  returned  to  the  Option  mode.  If  the 
user  answers  ”0"  to  the  question  about  using  the  default  values 
for  the  settling  time  calculations,  then  this  example  would 
continue : 

DO  YOU  WISH  TO  USE  DEFAULT  VALUE  OF  WITHIN  2Z 

OF  THIS  VALUE  FOR  CALCULATION  OF  SETTLING  TIME? 

(1  FOR  YES/0  FOR  NO)>0 

ENTER  UPPER  BOUND  FOR  SETTLING  VALUE 

>40 

ENTER  LOWER  BOUND  FOR  SETTLING  VALUE 

>39 

SETTLING  TIME  FOR  Yl-4 . 399999999999 
If  the  response  does  not  settle  down  within  either  the  user 
specified  band,  or  within  the  +2Z  of  final  value  range,  then  the 
message 

SIMULATION  DID  NOT  REACH  A  SETTLING  TIME 
is  displayed,  and  the  user  is  returned  to  Option  query. 

Almost  all  of  the  added  code  for  Option  #28  is  contained  in 
the  last  overlay  in  MULTI,  OVERLAY (2  1 , 0 ) . 

Suggested  Addition  To  Option  #28 

This  section  describes  program  code  that  would  calculate  and 


display  the  maximum  positive  aad  negative  actuator  rates. 


First,  f'-'o  common  blocks  have  to  be  added  to  OVERLA  Y(  2 1 , 0 )  : 
COMMON  /B  6  /  DMATR IX , ACT , SEN 
COMMON  /B  13A/  UP(101,11) 

Next,  right  before  the  common  block  declaration  area  add: 

CHARACTER  DMATR IX, ACT, SEN 
At  the  end  of  the  DIMENSION  statements  add: 

DIMENSION  RATEMAX( 10) ,RATEMIN( 10) 

At  the  end  of  the  DATA  statements  add: 

DATA  RATEMAX,RATEMIN/20*0 . 0/ 

Add  the  following  to  the  indicated  DO  loop: 

DO  1400  I-l.P 

DO  1400  J-l.LCOUNT 

IF(J.NE.l)  THEN 

RATE*  (  DP  (J  ,1+1  )~UP( J-l  ,I+1))/(UP(J,1  )-UP(  J-l ,  1 )  ) 

IF (RATE . GT.RATEMAX( I ) )  RATEMAX( I ) -RATE 
IF (RATE. LT.RATEMIN(I ) )  RATEMIN( I )-RATE 

END  IF 

1400  CONTINUE 

RATEMAX  and  RATEMIN  are  used  to  store  the  maximum  positive  and 
maximum  negative  values  of  the  rates,  respectively.  RATE  is  used 
to  calculate  each  rate,  with  UP(J,I+1)  storing  the  current  value 
of  the  actuator  deflection  and  UP( J-l , 1+1 )  storing  the  previous 
value  of  actuator  deflection.  UP(J,1)  and  UP(J-1,1)  store  the 
time  of  the  current  and  the  time  of  the  previous  actuator 
deflections  respectively.  To  display  the  results,  the  following 


code  Is  added  Co  Che  end  of  Che  overlay: 

DO  1600  l-l,P 

IF(RATEMIN(I) .LT.O)  THEN 
IF(ACT.EQ.'Y')  THEN 

PRINT*, 'MAXIMUM  NEGATIVE  ACTUATOR  RATE  # ' , I , '-' ,RATEMIN( I ) 
ELSE 

PRINT*, 'MAXIMUM  NEGATIVE  CONTROL  INPUT  RATE 
RATEMIN(I) 

ENDIF 

ENDIF 

IF(ACT.EQ.'Y')  THEN 

PRINT*, 'MAXIMUM  POSITIVE  ACTUATOR  RATE  # ' , I , ' ,RATEMAX( I ) 
ELSE 

PRINT*, 'MAXIMUM  POSITIVE  CONTROL  INPUT  RATE  #', I,'-' 
RATEMAX(I) 

ENDIF 

PRINT*,  '  ' 

1600  CONTINUE 

This  code  correcCly  informs  Che  user  Chac  Che  races  are  for 
conCrol  inputs  when  Chere  are  no  actuators. 

This  code  has  been  tested  on  a  copy  of  MULTI  and  seems  to 
work  correctly. 


Suggested  Addition,  Option  #38 

This  section  lists  the  code  that  would  implement  Option  #38, 
an  option  to  toggle  between  aodified  spline  and  linear  curve 
fitting  for  the  CALCOMP  plots.  For  aore  information  on  how  this 
code  works,  the  user  is  referred  to  the  CALCOMP  USER'S  GUIDE 
available  froa  either  the  AFIT  Engineering  Building  (Bldg  640)  or 


from  the  ASD  Computer  Center,  Building  676. 

The  first  step  is  to  add  the  following  to  the  common  block 
statements  in  the  main  section  of  the  program: 

COMMON  / B  13A/  UP( 1 01 , 1 1 ) , LINEFIT 
and  also  add  the  following  to  the  DATA  list  in  the  same  section: 
DATA  LINEF IT/-1 / 

Next,  add  the  following  to  0VERLAY( 1 0, 0) : 

First,  to  the  common  block  list: 

COMMON  /B  13A/  UP(101 , 11 ) .LINEFIT 
To  the  listing  of  options: 

PRINT*,  '  38.  TOGGLE  CALCOMP  PLOT  CURVE  FITTING" 

To  the  sequential  option  code  listings 
C  OPTION  #38 
2038  LINEFIT-LINEFIT*(-1) 

PRINT*, "CALCOMP  PLOT  CURVE  FITTING  IS' 

IF (LINEFIT. EQ.l)  PR INT* , 'L INEAR" 

IF (LINEFIT . EQ  •  *"1  )  PR INT* , 'MODIFIED  SPLINE" 

NOPT-O 
GO  TO  8010 

Then  add  the  following  to  OVERLA Y( 1 2 , 0) : 

To  the  common  block  list: 

COMMON  /B  13A/  UP ( 1 01 , 1 1 ), LINEFIT 
In  two  places  in  the  code  section  replace 
CALL  FLINE  (XAXI S , YAXIS , -N , 1 , 1 , 72 ) 

with 

CALL  FLINE  (XAXIS, YAXIS , N*LINEFIT, 1 , 1 , 72) 

This  option  has  been  added  to  a  copy  of  MULTI  and  tested 
with  no  problems.  Option  #38  acts  like  a  switch,  and  is  designed 


above  ) 


to  be  used  whenever  the  modified  Option  #28  (with  rates, 
displays  rates  that  the  user  thinks  are  too  high  for  the  modified 
spline  fit.  These  two  options  were  used  in  the  author's  thesis 
to  prevent  actuator  deflection  plots  from  assuming  non-functional 
shapes. 


APPENDIX  D 

SYSTEM  MATRICES,  CONTROLLER  MATRICES,  AND  DESIGN  PARAMETERS 
FOR  0.6  MACH,  SEA  LEVEL,  AND  FOR  0.9  MACH,  30,000  FEET 


Introduction 

This  appendix  lists  the  system  matrices,  the  design 
matrices,  and  the  design  parameters  for  the  two  flight  conditions 
not  discussed  In  detail  In  the  text  of  this  thesis,  0.6  mach  at 
sea  level,  and  0.9  mach  at  30,000  feet.  The  complete  system 
matrices  In  the  form  of  equation  (48)  from  Chapter  III  can  be 
found  In  Jon  Bauschllcher's  thesis  (Ref  2). 
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